DD6  FILE  COPY, 


* y 


^W«tt I 


WmgmmM 


$Bs&m 

lisps 


AFOSR-TR-  7 9-0960 

SYSTEMS  TECHNOLOGY,  INC 


J?*4  SO  NAWTNOINf  lOUltVAID  • HAWTKOtNI,  CALIFORNIA  *0330 


Final  Scientific  Rcoort  No 


MODELING  BIODYUAMIC  EFFECTS  OF  VIBRATION 


Fifth  Year  — Final  Scientific  Reoort 


August  1979 

Approved  fer  publio  r*lwu#i 


(B 


SYSTEMS  TECHNOLOGY.  INC. 

s°l  Til  II  WVT1IORM-;  M)«  I.KVVRI)  • 1 1 WVTHORNE . CALIFORNIA  H02S0  * PH<X\E  12131  A7D-22H I 

\*f  i * <■ 


Final  Scientific  Report  No. 


‘ MODELING  BIODYNAMIC  EFFECTS  OF  VIBRATION 
Fifth  Year  — Final  Scientific  Report 


H.  R.  Jex  • 
R.  E.  Magdaleno 


r-r-vriFrc  RsisiF.cn  (a?sc.j, 


clinical  Information  Officer 


f:.d  IS 


August  1979 


Approved  for  rj 

v ,r.  ri  •.-.tion  unlinitett 


Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research/ 

AFSC,  United  States  Air  Force,  tinder  Contract  No.  FUU620-73-C-0075 . 

The  United  States  Government  is  authorized  tc  reproduce  and  distri- 
bute reprints  for  governmental  piurposes  notwithstanding  any  copy- 
right notation  hereon. 

79  09  14  077 

lilt. .A.  II  HI- 1 It  I.  Nit’ll  * N't*  \ | N VIEW.  CU.IMHIMA 


SECURITY  CL  ASSl  PIC  ATION  of  THIS  PAGE  Oats  Entered) 

REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTION'S 
BEFORE  COMPLETING  FORM 


REPORT  NUMBER 


' I« — IlUL6.C*l4JS«A^Wa_ 


!2.  GOVT  ACCESSION  NC.j  3.  RECIPIENT'S  CATALOG  HUMBER 


fv«T>f-wi»osTi  ptwct  ':o»*Af-3 


>T*ariTT  i hBl 


l—1 l— (Fifth  Year#  — » r*ml'  io  twpoyt) 


»f:-1^37-5  J£ 


3.  R E P ORT  f,ruW8iJL/ 


■M — 

KJ  R.  fizx  ■ ) \ 

r ;•/«■  E. /Magdalene/ 

. "9.  PERFORMING  ORGANIZATION  NAME  ANO  ACORESS 

Systems  Technology,  Inc/ 

1 3766  South  Hawthorne  Boulevard  Qu) 

Hawthorne,  California 2Q,25fl Z^rz 

U.  CONTROLLING  OEFtCE  NAME  AND  ADDRESS  ■/ 

Force  Cfv.ce  o(  Sde^Wf.c  Qc^arcL_(/Vt/  {// 
Art=b 

— u^^Tnamr.  ^EitCY  AM£  4 ADDRESS (if  dlftoront  from  Controlling  Otficm) 

16-  DISTRIBUTION  STATEMENT  (ot  fhi»  Report) 


i,  CQMTfiM>€-*"»4*"GRAN7  N U M 8 E R/ 1 ) 

u i?>C-C>cn  5V 

ibT  PROGRAM  EL  EMC  NT.  PROJECT.  TASK 
AREA  6 WORK  UNIT  NUMBERS 

to\\C&*cZ?l 
«33>~  1 o)  }/  A SL.  ( 

u RfctNartT  a*»g — 

Augwt  *979/ 

li  NUMBER  OP  PAGES 

H£i_ 

IS.  SECURITY  CLASS,  (ot  thla  report) 


1$«.  DECLASSIFICATION-  DOWNGRADING 
SCHEDULE 


Approved  for  public  rolcage;  diofrrib 

OD/irtefi 


y die%ribnt ion  unlimited. 


•'?.  DISTRIBUTION  STATEMENT  to)  th 


/ v ^-7?-///^, 


i*.  supplementary  notes 


M9.  AEY  ROROS  rContlmrt  on  etde  tl  neceeeery  and  Identity  by  block  number) 


Vibration  Visual  acuity,  dynamic 

Human-operator  BIODYN-78  Computer.  Program 

Transmissability 
Math-models 

Vestibulo-ocular  reflex 

?vK>^A0STR  ACT  fContimro  on  %t4o  H n«c«i«ary  antf  Identity  by  block  mm »6#r> 

The  biomechanical  feedthrough  of  vertical  and/or  fore-aft  vibration 
from  a seat,  to  (and  through)  various  parts  of  the  torso,  limbs,  hands,  and 
head,  .to  a spring-restrained  control  stick  has  been  modeled  over  a 5-year 
effort  using  prior  Air  Force  and  other  data  bases.  This  final  report 
summarizes  the  project  results  and  provides  a User's  Guide  to  BIODYN-78, 
sm  interactive,  remote-access,  digital  computer  program  for  exercising  the 
model. ~ — - 


DD  1 JAN*!  .1473  COITION  O#  I NOV  ••  IS  obsolete 


MS' 


Unclassified  / 

** CURITY  CLASSIFICATION  Of  THIS  RAGE  (When  DeufEr^Cd) 


SECU  PlTV  CLASSIFICATION  OF  THIS  P AGEfHTi./.  On(«  Enfrmd) 


■The  model  involves  77  postural  and  physiological  parameters  to  define  a 
broad  range  of  seating/control/display  arrangements  from  supine  to  erect.  An 
internal  routine  automatically  trims  the  Torso/Head/Limb  system  at  the  desired 
angles  and  computes  linearized  equations  of  motion,  which  are  solved  for  their 
eigenvalues  and  resulting  transfer  functions  using  highly  refined  computing 
routines.  In  distinction  from  most  past  models,  this  model  represents  several 
degrees-of- freedom  simultaneously,  e.g.,  vertical  and  fore/aft  shoulder  motion, 
angular  head-bobbing,  constrained  elbow  motion,  stick-grip- interface  force, 
muscle  forces,  and  stick  angular  motion,  among  many  others . These  features 
make  the  BIODYN-78  model  well  suited  for:  routine  computations  of  biomechanical 
transmissibilities,  (e.g.,  those  required  as  inputs  to  other!  programs  (PIVIB, 
USAM)  for  pilot/vehicle  closed  loop  performance  calculations^) , matching  of 
test  data-,  systematic  analysis  of  variations  in  body  parameters  (e.g.,  to 
explore  population  extremes),  optimization  of  control  stick  "feel"  parameters 
(i.e.,  to  analyze  Pilot  Induced  Oscillations),  or  to  compute '.control  transients 
due  'to  accelerations  imposed  on  the  pilot  by  maneuvers.  j 

The  fifth-year  effort  covered  two  objectives.  The  first  objective  of  the 
FY  1978  program  was  to  evolve  a model  for  the  effects  of  head  vertical  plane 
vibrations  on  eye  and  image  motions;  (as  influenced  by  the  vestibule-ocular 
and  fixation  reflexes).  This  was  met  by  the  model  detailed  in  the  Fourth 
Interim  Scientific  Report  and  summarized  in  Appendix  A*  herein. 

The  second  objective  was  preparation  of  the  User's  Guide  (included  here  as 
Appendix  Bf),  and  the  formal  demonstration  of  BIODYN-78  operation  using  a 
remote-access  terminal  in  the  Aerospace  Medical  Research  Laboratory  (at  WPAFB, 
Ohio). 


* Progress  in  Measuring  and  Modeling  the  Effects  of  Low  Frequency  Vibration 
■ on  Performance,  Paper  No. .232 


t "User's  Guid-?  to  BIODYN-78,  An  Interactive  Computer  Program  for  Modeling  • 
Biodynamic  Feedthrough  to  a Pilot's  Hands  and  Head,  TR-1 037-2 
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A.  BACKGROUND 


1 . Scope 

This  final  (fifth-year)  report  on  the  project  briefly  reviews  the 
objectives  and  work  covered  in  the  prior  four  and  one-half  years,,  and 
describes  in  more  detail  the  work  done  in  the  last  half  of  the  fifth  year, 
i.e.,  May-December  1978.  (The  last  Interim  Scientific  Report  covered 
work  into  half  of  FY>  ' 78  because  that  offered  a more  logical  reporting 
break  point. ) The  main  text  describes  the  results  and  recommendations, 
while  two  appendices  give  more  details: 

• Appendix  A — Contains  a status  report  on  the  state-of- 
the-art  in  modeling  biodynamic  effects  on  aircrew  manual- 
task  performance  (Ref.  l).  It  was  an  invited  paper 
presented  at  the  AGARD  Aerospace  Medical  Panel  Specialists 
Meeting,  Paris,  Nov.  1973. 

• Appendix  B — Includes  a brief  "User’s  Guide  to  BIODYN-78, 

An  Interactive  Computer  Program  for  Modeling  Biodynamic 
Feedthrough  to  a Pilot's  Hands  and  Head”  (Ref.  2). 


2.  Objectives 

The  various  phases  of  this  research  program  have  the  primary  objective 
of  developing  analytical  dynamic  models  of  the  biomechanical  effects  of 
vibration  on  a pilot's  control  behavior  and  performance,  which  could  be 
used  by  aerospace  system  engineers' and  applied  researchers  in  this  field. 
The  prior  criteria  set  by  APC3R  for  such  engineering  models  are  that  they 


should: 


a.  Be  based  on  a sound  empirical  data  base. 

b.  Capture  the  essential  phenomena  of  interest  to  the 
users,  without  excessive  restrictions  or  details. 

c.  Be  compatible  with  more  detailed  "microscopic"  model 
results. 

d.  Be  put  in  a forn.  usable  by  others j at,  above,  and 
below  in  the  hierarchy  of  applications. 
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The  prior  four  years'  work  attempted  to  meet  the  first  three  objectives, 
and  this  final  phase  meets  the  last  objective. 

The  specific  work  items  of  the  fifth  year  contract  were  to: 

• Add  a detailed  model  for  the  effects  of  vertical 
vibration  on  the  head-eye  motions  (including  both 
’'fixation-reflex"  loops  and  "vestibulo-ocular  re- 
flex” loops). 

• Put  the  1978  biomechanical  model  (prior  to  eye- 
motion  effects)  in  a time-shared  computer  program 
accessible  to  any  inheres  ed  and  qualified  user, 
including  output  of  biomechanical  transfer  func- 
tions r quired  as  inputs  to  the  ARML/BEN  "PIVIB-76', 
tracking  performance  program. 

• (Carried  over  from  FY  *77)  on  an  informal  basis, 
assist  AMRL/BB  with  the  analysis  of  the  beadrbobbing 
experiments  of  Dr.  P.  Werner  for  possible  use  in 
validating  the  BIODYN  model  at  low  frequencies. 

Before  describing  the  results  of  this  year's  work,  we  will  present  a 
brief  overview  of  the  approach  used  and  the  work  performed  during  various 
phases  of  this  project,  to  give  a proper  perspective  to  this  final  report. 

3 . Approach 

Our  approach  is  to  use  control  systems  theory  and  dynamic  analysis  1 
techniques  to  assist  in  the  measurement  and  modeling  of  pilot  behavior  (i.e., 
frequency  responses,  describing  functions,  dynamic  feedthrough  functions, 
etc.)  and  performance  (e.g.,  rms  errors,  coherency  to  forcing  functions, 
stability  margins,  visual  impairment,  etc.),  under  various  types  of  vibra- 
tion (sinusoidal,  random,  vertical,  lateral,  etc.).  Hie  empirical  data  base 
relies  on  prior  published  experimental  data  and  on  new  experimental  data 
provided  by  the  6730  AMRL/BB  Branch  at  Wright-Patterson  Air  Force  Base 
(under  the  overall  direction  of  Dr.  H.  von  Gierke). 

The  biomechanical  models . derived  use  lumped-parameter  approximations 
to  the  dominant  body/limb/head  configurations  likely  to  be  encountered  in 
pilots  of  aircraft,  where  an  erect  or  semi-supine  posture  is  present,  often 
with  some  redundant  constraints  on  motions,  due,  to  seat  restraints,  arm- 
rests, or  control-stick  gripping.  The  motions  have  been  separated  into  a 
SR- 1037-1  2 


lateral  set  (sway  inputs)  and  a vertical-plane  set  (heave  and  surge  in- 
puts). Only  the  vertical-plane  set  has  been  carried  beyond  the  Dasic 
modeling  stage,  because  of  its  dominant  importance  for  high-load- factor 
maneuvers  and  because  of  scarcity  of  data  for  the  lateral  set  (Ref.  23). 

Our  models  differ  from  the  three  most  commonly  used  types  of  bio- 
mechanical models,  which  include: 

a.  One-input,  one-direction-output  models  such  as  the 
vertical  input,  vertical  output  whole-body  vibration 
models  used  decades  ago  by  Coermann  (Ref..  3)  expanded 
by  von  Gierke  and  Goldman  (e.g..  Ref.  k)  and  Payne 
(Ref.  5)  and  most  recently  extended  to  nonlinear 
effects  of  high  Gz  by  Mertens  (Ref.  6).  These  models 
are  only  intended  to  represent  specific  outputs  in 
the  same  direction  as  the  input  (i.e.,  vertical  shoul- 
der or  head  motion  outputs),  and  they  cannot,  in  gen- 
eral, simultaneously  represent  multiraple  outputs  (e.g., 
head  angle,  arm  angle,  shoulder  heaving,  etc.)  or  mul- 
tiple inputs  (e.g.,  heave  and  surge). 

b.  Open-chain,  planar,  passive  models  such  as  the  numerous 
body-trajectory  models  developed  for  crash  research  at 
Calspan  and  elsewhere  (Refs.  7 and  8)  and  applied  by 
Kaleps  (Ref.  9)-  3uch  models  and  the  efficient  canonicax- 
form  equations  used  therein  are  excellent  for  integrating 
large,  unimpeded  motions  but  are  not  appropriate  for  con- 
strained or  high-frequency  vibrational  motions. 

c.  Nonlinear^  continuous-media,  or  detailed- joint  models  such 
as  the  spinal-column  models  of  Ref.  10,  and  neck  models  of 
Ref.  11.  These  very  derailed  models  of  every  joint,  and 
nonlinear  visco-elastic-  connecting  elements  are  the  "micro- 
scopic" models  referred  to  earlier.  They  can  represent  all 
the  possible  modes,  degrees  of  freedom  and  inputs,  but  are 
too  expensive  to  exercise  for  most  routine  problems  of 
pilot-vehcile  optimization,  where  only  the  dominant  modes 
need  be  modeled. 

The  present  models  are  aimed  at  pilot-vehicle  control  interface  problems 
in  biodynamic  environments  (e.g.,  tracking  an  enemy  fighter  at  high  load- 
factor,  wherein  transonic  wing-buffeting  is  shaking  the  aircraft  randomly 
in  the  1-20  Hz  frequency  range).  To  meet  these  needs  our  models  have  che 
following  features,  many  of  which  differ  from  those  preceding:, 
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a.  "Homologous  elements,"  capable  of  representing  simultaneously 

the  several  degrees  of  ..’“eedoni  of,  the  dominant  modes  for 
the:  torso/arm/hand/s  .ick  and  torso/neck/head/ eye/ display 

chains.  These  chtii..;  involve  closed  kinematic  and  neuro- 
muscular loops,  and  their  models  can  accept  inputs  in 
several  degrees  of  freedom,  simultaneously. 

b.  All-posture,  nonlinear  kinematics  (supine  to  seated); 
linearized  about  a given  equilibrium  posture  (appropriate 
for  the  relatively  small  vibrational  motions  present). 

c.  A physical  equations-of-motion  approach,  in  ■which  a given 
(lumped)  parameter  appears  only  once  in  each  equation. 

This  greatly  facilitates  parameter  matching  when  fitting 
a model  to  empirical  data. 

d.  Element  masses  are  retained  (not  divided  into  the  damping 
and  compliance  as  is  customary  in  modal  equations),  so 
that  some  mass  elements  can  be  finite  or  zero,  without 
blowing  up  the  equations. 

e.  Extensive  use  of  linearized  transfer-function  formats 
and  factored  roots  is  used  to  present  the  resulting 
dynamics  in  a form  most  useful  to  systems  dynamicists 
(e.g.,  as  frequencies  and  damping-ratios  of  various  trans- 
fer function  numerator  and  denominator  factors). 

f.  Provision  of  active  adjustable,  neuromuscular  loops 
for  the  limb,  including  "filtered"  spindle  and  golgi 
feedback  effects.  These  influence  the  limb's  apparent 

1 impedance  to  vibrational  relative  motions  between 

torso  and  hand  or. head. 

Further  details  of  our  approach  are  given  in  Appendix  A (i.e..  Ref.  l).  Our 
approach  to  supplying  the  numerous  coefficients  involved  is  to  use  actual 
dimensions,  masses,  and  rotary  inertias  whenever  possible,  based  on  cadaver 
data  or,  in  some  cases,  on  direct  measurement  (e.g.  f forearm  weights  for 
light  and  heavy  test  subjects).  The  compliance  and  damping  coefficients 
are  then  selected  by  a,  combination  of  prior  data  and  model  fitting  to 
match  particular  modes  of  empirical  data. 

To  fit  such  complex  models  (over  20th  order)  to  various  experimental 
results,  a special  technique  was  developed,  the  "Multiloop  Model  Fitting 
Procedure,"  which  is  described  in  our  Second  Interim.  Scientific  Report.  . 

(Ref.  15).  Suffice  it  to  say  here  that  this  technique  allows  the  matrix 
expansion  of  a given  transfer  function  numerator  and  denominator,  in  ' \, 
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terms  of  summed  matrix-products  with  each  parameter  which  is  to  be 
varied,  thus  permitting  a far  more  efficient  parameter  optimization 
procedure.  The  frequency  response  vectors  at  fifteen  frequencies  in 
three  different  degrees-of-freedom  are  simultaneously  fitted  by  a 
weighted-least-squares-criterion.  Over  the  program,  roughly  one-quarter 
of  our  effort  was  spent  in  developing  and  using  these  model  fitting 
procedure's  to  refine  and  validate  the  models. 

4.  Chronology  and  Reports  Issued 

Figure  1 shows  the  chronology  of  this  APOSR  program,  in  relation  to 
the  related  experimental  work  at  the  USAF  Aerospace  Medical  Research 
Laboratory,  Biodynamics  Branch,  at  WPAFB  (under  Dr.  von  Gierke)  and  at  the 
Office  of  Naval  Research/Human  Factors,  Inc.,  Ship  Motion  Generator' 
Facility,  then  at  Goleta,  CA  (under  cognizance  of  Dr.  Arthur  Callahan 
of  ONR  and  the  USN  Surface  Effect  Ship  Project).  Table  1 summarizes 
some  of  the  reports  prepared  under  the  Air  Force  vibration  programs. 

The  first  two  years  of  this  program  covered  a detailed  reanalysis  of 
systematic  vertical  and  sway  vibration  experiments  using  sinusoidal  vibra- 
tions. Simplified  transmissibility  models  were  developed  for  torso-limb 
for  both  vertical  and  sway  inputs.  Performance  data  from  both  sinusoidal 
and  new  quasi-random  vertical  vibration  data  were  correlated  and  shown  to 
be  primarily  dependent  upon  the  "characteristic-frequency"  of  vibration 
(frequency  of  upward  axis  crossings  of  vertical  acceleration)  and  its 
amplitude  (Refs.  1^  and  15). 

The  third  year  covered  intensive  development  of  the  Multi loop  Fitting 
Procedure,  inclusion  of  a one-element  head  bobbing  made  to  the  vertical 
model  and  addition  of  an  arm-rest  constrai..;  on  limb  fore-aft  sliding  • 
motion.  Although  it  complicates  the  equations  the  last  iteto  is  important 
because  numerous  pilot- vehicle  situations  involve  an  armrest  (e.g.,  side- 
sticks  and  center-sticks  with  the  arm  rested  or.  the  knee) . Extensive  model 
fitting  to  the  previous  transmissibility  data  was  performed,'  and  needed 
refinements  were  revealed  (e.g.,  a better  head/neck  ncdel.  Ref.  16). 
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Figure  1.  Chronology  of  Related  STI  Projects  for  AMRL/BBand  AFOSR 
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The  fourth  year  included  developing  the  indicated  passive  neck/head 
model,  and  its  validation  by  some  experimental  data.  We  assisted  Dr.  Peter 
Werner  at  A.MRL/BB  in  planning  a series  of  head-bobbing  experiments  on  the 
WPAFB  Vertical  Accelerator  facility  to  cover  the  scarce-data  range  of 
0.5  to  5 Kz  and  O.i  to  0.3  Gz.  Meanwhile,  some  head-point-of- regard  data 
taken  as  a byproduct  of  a U.S.  Navy  program  (on  th_  "habitability"  of 
high-speed  Surface  Effect  Ships  bouncing  over  rough  waters)  became  avail- 
able. Through  the  cooperative  efforts  of  AFOSR,  ONE,  arid  the  David  Taylor 
Naval  Ship  Research  and  Development  Center  (DTNfRDC),  several  digital  tapes 
of  head  notion  under  various  random  input  and  postural  conditions  (resting, 
reading,  tracking)  were  processed,  providing  a rich  and  reliable  data  base 
which  has  not  yet  been  fully  analyzed.  The  refined  model  will  fit  head- 
point-of-regard  motions  in  the  ''passive"  and  "reading"  cases  but,  as  noted 
in  last  year's  Interim  Scientific  Report  (Fig.  3 of  Ref.  19) , there  is 
evidence  for  active  neuromuscular  head  control  in  the  "head-tracking"  pos- 
tural condition.  This  remains  for  future  wcrk.  During  this  fourth  period 
we  extended  the  equations  to  cover  serisupine  postures  and.  axial  (fore-aft) 
acceleration  inputs . 

It  became  apparent  during  this  time  that  visual  performance  decrements 
due  to  vibration  were  dependent  on  the  eye's  image  motions  induced  by  head 
motions,  and  that  the  eye's  sopni3ticated  control  system  would  require  an 
equally  sophisticated  model  to  ccpe  with  the  recently  established  motion 
properties  of  the  "fixa* ion-reflex"  (more  properly,  a closed  target-tracking 
loop)  .and  the  "vestibulo-ocular  reflex"  (inertial-motion  compensation  loop). 
Thi3  model  «*as  begun  in  the  fourth  year  and  finished  during  the  first  por- 
tion of  th®  fifth  year,  so  it  was  reported  in  last  year's  ISR.  The  eye- in- 
head model  was  developed  for  vertical  angular  motions  of  the  head,  tout  is 
also  applicable  to  head  yawing  motions  with  suitably  adjusted  parameters. 

It  seems  to  correlate  well  with  some  recent  English,  data  on  visual  acuity 
under  vertical  vibration,  and  offers  a final  link  in  the  chain  between  bio- 
dy..-4flic  inputs  to  the  body  and  decrements  in  visual  performance.  Further 
experimental  data  and  pursuit  of  the  implications  of  this  model  on  the  use 
of  helmet-mounted-display3  seems  warranted.  (See  the  discussion  in  last 

year's  ISR,  Ref.  19,  and  ir.  Appendix  A herein). 

This  summary  brings  us  to  the  final  half  of  the  fifth  and  last  year, 
uhich  is  reported  on  herein. 
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B.  WORK.  DURING  FY  1978 

1 . Completion  of  Rbdel  for  Vibration  Effects  on 
the  Eye  Control  System 

As  notea  above,  the  eye-control  model  begun  in  FY  1977  was  finished  dur 
ing  the  first  half  of  FY  1978  and  was  included  as  Appendix  B of  last  year''*, 
summary  report  (Ref.  l8,  included  in  Ref.  19).  An  overview  of  that  model 
and  some  key  results  Is  given  herein  in  Appendix  A (Ref.  l).  Consequently, 
we  will  not  describe  It  here.  It  is  being  submitted  to  an  archival  jour- 
nal. 


The  state-of-art-  survey  paper  included  here  as  Appendix  A was  an 
invited  presentation  at  the  AGARP  Aerospace  Medical  Panel  Specialist's 
Meeting  in  Paris  during  November  1979-  At  the  suggestion  of  Dr.  von  Gierke 
of  AMRL/BB,  it  describes  the  approach  and  typical  results  for  work  on 
vibration  effects  done  both  by  Systems  Technology,  Inc.  (predominantly  on 
biomechanical  and  visual  performance  modeling)  and  by  Bolt,  Beranek  and 
Newman  (predominantly  on  operator  tracking  performance). 

During  visits  to  some  vibration  laboratories  in  Europe  subsequent  to 
that  presentation*  the  relevance  of  our  biodynamic  models  to  a number  of 
AGARD-related  experiments  and  problems  was  discussed: 

• At  DFVLR  - Institute  for  Aviation  Medicine  (Bad  Godesburg, 

West  Germany),  prs.  H.  L.  Vogt  and  H.  Martens  were  interested 
in  reconciling  their  uni-directional  vertical  vibration  model 
with  the  STI  model.  The  appropriate  combinations  of  STI  and 
DFVLR  parameters  was  discussed,  and  a further  conamrication 
has  been  received  from  Dr.  Mertens.  If  we  can  establish  the 
proper  relationships  between  these  two  models  and  parameters, 
then  the  STI  model  can  use  the  extensive  data  of  DFVLR  on  the 
covariation  of  various  torso  arid  limb  compliances  with  steady 
Gz  level  to  extend  our  model  to  simulate  the  rapid-onset  of 
load  factor.  Such  experiments  are  currently  of  great  interest 
at  AMRL/BB. 

• Max  Planck  Institute  for  Agricultural  Research  (Bad  Kreuznachj 
We3t  Germany).  Dr.  H.  Dupuis,  who  has  had  a significant  role  ■ 
in  establishing  the  ISO-Ride  Quality  Criteria,  is  doing  work 


* These  visits  were  not  made  on  APOSR  funds. 


on  the  effects  of  the  NATO  Main  Battle  Tank  vibrations  on  the 
gunner's  perceptual  problems.  Hie  gunner  holds  his  head 
against  a hairi-padded  rest  to  minimise  whole-body  vibration 
effects*  but  higher  frequencies  cause  the  eyeball  to  have 
rotary  resonances  with  resulting  image  blur.  Dr.  Dupuis 
noted  that  the  eye-control-resonance  frequency  and  damping 
used  in  our  eye-control  model  (estimated  from  earlier  labora- 
tory  data  on  eye  saccades)  coincided  closely  with  his  recent 
experimental  findings  using  a corneal-reflective  TV  tracking 
scheme.  He  was  also  interested  in  our  posturally  complete 
model  for  representing  several  degrees-of- freedom  simultaneously, 
as  an  improvement  over  the  one-direction  models  commonly  used, 
and.  he  supplied  a number  of  useful  references  on  limb  and  hand 
vibration  effects. 

9 Krupp  Hospital's  Neurological  Clinic  (Essen,  West  Germany). 

Dr.  Thomas  H.  Brandt  is  doing  work  for  the  NASA  Space  Shuttle 
Vestibular  Experiment  and  in  clinical  evaluation  of  vestibular 
system  deficiency.  They  were  contenqiLating  a device  for  creat- 
ing sudden  angular  impulses  to  the  head  to  check  the  operation 
of  the  ve st ibulo -ocular  reflex  before  the  fixation  reflex  could 
re-stabilize  the  line  of  sight.  They  had  been  having  problems, 
which  our  model  explained.  Discussion  of  the  model  implications 
suggested  an  alternative  angular  vibration  approach,  with  inputs 
in  the  frequency  rarge  which  exhibits  low  sensitivity  to  fixation 
loop  dynamics  ,but  high  sensitivity  to  the  vestibular  loop  cross- 
feeds. Development  of  this  clinical  application  is  recommended. 

• RAF  Institute  of  Aviation  Medicine  (Pamborough,  England). 

Dr.  /Vi an  Benson '3  colleague  Dr.  Graham  Barnes  has  done  much 
work  on  measuring  and  modeling  the  vestibulo-ocular-reflex 
and  nystagmus  effects  of  rotary  head  motions  (our  model  builds 
3trongly  on  this  work).  He  has  offered  to  review  Ref.  18, 
and  to  supply  3ome  more  recent  data  on  the  effects  of  visual 
3cene  complexity  cn  the  fixation  loop  gain,  specifically  to 
help  fill  in  a data-gap  mentioned  in  our  writeup. 

At  the  separate  RAE  Human  Engineering  Division,  under  the 
direction  of  Mr.  Geoff  Allen,  his  colleague,  il.  F.  Rowland 
agreed  to  supply  some  essential  data  on  pilot  seat-posture 
and  dimensions  missing  from  his  detailed  studies  of  head 
vertical  motions  as  influenced  by  resting  posture.  Especially 
Interesting  are  the  "back-off- seat'*  Posture  which  produced 
dramatic  reductions  in  vertical  head  motions  under  vertical 
vibrations.  If  our  model  can  "explain"  these  independent 
data,  then  a powerful  validation  and  impetus  to  new  applica- 
tions will  have  been  obtained.  However,  as  of  th<  «nd  of 
thi3  contract  (December  1978)  the  information  had  not  yet 
arrived. 
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These  researchers  would  be  interested  in  using  the  BIODYN  racdel  to  help 
fit  and  codify  their  data,  if  it  were  available  in  a readily  exercisable 
form.  If  this  could  be  done  cn  an  information-exchange  basis,  then  it  would 
be  to  USAF's  advantage  to  work  out  some  arrangement  for  a time- shared-network 
version,  by  which  U3AF  would  gain  the  fitted  data  as  part  of  a much  broader, 
"fitted"  data  base  at  relatively  little  cost,  while  the  users  would  gain  the 
program's  advantage  at  no  software  development  cost.  We  will  return  to  this 
suggestion,  later. 

2.  Computerized  Transmissibility  Model  — 

"BIODYN- 1978” 

a.  Objectives 

As  noted  in  the  background  section,  one  criterion  of  a mature  math 
model  is  that  it  is  in  a form  readily  understandable  and  usable  on  the 
problems  it  is  intended  to  solve.  Some  of  the  intended  roles  for  this 
biomechanical  interface  model  are  as  follows: 

^ • To  guide  experimenters  on  the  most  useful  locations  and 

frequencies  for  vibration  measurements,  and  to  fit  data 
from  past  and  current  experiments. 

9 In  a case-by-case  "user-interactive"  mode  by  R and  D 
lab  personnel  to  develop  a part icular  model  posture 
or  parameter  ensemble,  or  in  a "batch-mode"  to  syste,- 
matically  map  out  a range  of  vibration  effects. 

• To  be  routinely  used  by  aerospace  development  engi- 
neers in  solving  practical  design  problems  such  as 
Pilot  Induced  Oscillations  (often  due  to  limb-stick 
"bobweight"  effects);  posture  effects  from  High- 
Acceleration-Cockpit  (HAC)  design  alternatives  (e.g. 
semisupine  seating,  head  or  torso  support  compliance); 
stick  location  and  arm-re3t  effects,  etc. 

• ' To  permit  flight  control  system  engineers  to  optimize 

combined  vehicle/aircrew  ride-qualities  and  visual 
performance  effects,  through  use  of  vibration  filters, 
dampers,  and  suppressors  (a  la  B-l);  or  for  display 
blur  minimization  through  choice  of  viewing  distance, 
angle  or  filtering,  etc. 

I 
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• To  provide  the  biomechanical  seat-to-hand  and  seat-to- 
head  transmissibility  and  limb/stick  impedance  transfer 
functions  required  by  the  AMRL/BBN  Pilot  VIBration 
(PIVIB)  tracking  performance  program  (Ref.  20).  Such 
biomechanical  functions  must  be  supplied  by  the  user  of 
PIVIB,  and  are  a natural  application  of  our  present 
model. 

The  time- shared  computerized  program  we  have  put  together  to  accomplish 
these  objectives  is  called  "BIOBYN-78")  the  number  specifying  the  model  year 
or  "vintage"  of  the  model,  since  additions  are  contemplated  in  the  future. 

For  example,  we  plan  to  include  the  head-eye  model,  described  earlier,  in  the 
BIODYN-79  version. 

b.  Approach 

The  BIODYN-78  model  includes  a chain  of  interacting  parallel  and  serial 
second-order  (mass/spring/damper)  elements,  with  complex  neural  ar.d  force 
feedbacks.  Details  were  given  in  earlier  reports  and  in  Ref.  17,  and  the 
mechanical  arrangement  is  shown  in  Fig.  2 of  the  Appendix  B herein  (Ref.  2). 
Figure  2,  reproduced  from  Ref.  l,  summarizes  the  functional  blocks  of 
BIODYN-78  and  denotes  the  number  of  parameters  involved  in  each.  It  also 
shows  the  relationship  of  BIODYN-78,  to  the' PIVIB  and  possible  future  BIODYN 
versions. 

The  equations-of -motion  for  BIODYN-78  have  been  written  in  & "second 
order,  element"  matrix  form  with  generalized  (all  angle)  parameters.  Any 
particular  set  of  postural,  neuromus cular , stick  and  display  dimensions  are 
then  input  as  a file  to  the  program.  The  program  internally  sets  up  the 
general  equations  and  evaluates  the  linearized  coefficients  for .perturbations 
about  the  equilibrium  posture.  The  program  then  evaluates  the  desired  trans- 
fer functions  (from  a "CHOICES " file)  via  Cramer's  rule  yielding  a demon - 
inator  and  various  numerators  of  rather  high  order  — often  reaching  20th 
order. 

Fairly  sophisticated  factorization  routines,  evolved  during  two  dec- 
ades of  fact< ring  aircraft/control  system  dynamics  equations,  are  used 
to  give  the  predominantly  second-order  poles  and  zeroes,  many  of  which 
are  lightly  damped.  Although  these  routines  are  considered  proprietary 
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STI  developments,  they  have  been  incorporated  into  the  accessible  BIODYN-78 
program  with  minimal  documentation. 

The  primary  outputs  of  BIODYN-78  are  the  factored  transfer  functions 
between  the  selected  input  (vertical  or  fore-aft  platform  acceleration  or 
stick  force)  and  any  of  the  36-odd  equation-of-motion  variables,  such  as 
head -bobbing  angle,  stick  motion  (control  transmissibility),  shoulder 
motion,  muscle  force,  grip  force,  etc. 

To  maximize  the  ease  of  accessing  and  using  BIODYN-78  we  have  put  at 
into  a commercial  remote-time-shared  computer  network  (Tymshare  Inc.), 
which  is  readily  accessible  anywhere  in  the  USA  and  in  certain  parts  of 
Europe  (via  satellite). 

The  use  of  Tymshare  for  BIODYN-78  has  several  advantages: 

• It  makes  BIODYN-78  potentially  available  to  any  user 
with  a low  speed  terminal  ( TO  or  30  cps)  and  acoustic 
coupler,  who  has  or  can  get  a Tymshare-Inc./STI  operat- 
ing agreement.  Simplified  procedures  for  accessing 
BIODYN  are  being  worked  on  (write  the  authors  for 
updated  information) . 

• Tymshare 's  powerful  EDITOR  is  available  for  modifying 
existing  files  or  creating  new  Input  Files  of  parameters. 
Excellent  User 1 s manuals  are  available  for  using  Tymshare 
in  the  simplest  possible  ways,  such  that  a non-programmer 
can  easily  operate  the  program. 

• Costs  of  re-writing  and  documenting  our  extensive  Tym- 
share-oriented  software  were  minimized,  as  none  of 
this  is  necessary  to  use  BIODYN-78. 

Although  BIODYN-78  is  currently  resident  in  the  Tymshare  Inc.  network 
for  the  reasons  given  above,  it  has  recently  been  adapted  to,  and  test 
run  on  Control  Data  Corporation's  CDC  time-shared  network.  For  the  same 
inputs  identical  outputs  resulted. 

Running  times  are  comparable  for  Tymshares  PDP-HO  or  CDC's  CYBER  175» 
The  main  advantage  in  using  Tymshare  Inc.  as  the  residence  for  BIODYN-78 
is  the  ready  availability,  therein,  of  STI's  Unified  Servo  Analysis  pro- 
gram (Ref.  21)  for  manipulating  and  plotting  the  results. 


— *■ 


c . Overview  of  Procedures 

A simplified  flow  diagram  for  the  procedures  used  in  applying  BIODYN-78 
is  given  in  Fig.  3.  The  User  accesses  Tymshare,  Inc.  (with  suitable  pass- 
words, etc.)  and  asks  for  one  of  the  already  filed  sets  of  parameters  for 
various  postures  and  situations  as  listed  in  the  User's  Guide  Catalog;  or 
he  can  create  new  ones  from  scratch.  If  some  of  the  parameters  need  to  be 
changed,  he  uses  Tymshare's  EDITOR  to  change  the  appropriate  numbers  on 
given  lines  of  the  selected  INPUT  File  (and  renames  it  accordingly),  or  he 
can  use  the  selected  set  without  modification.  This  process  could  take 
from  one  minute  to  one  hour  depending  on  the  number  of  parameters  to  be 
changed  and  checked. 

The  user  then  identifies  the  desired  transfer  functions  via  a "CHOICES" 
file  wherein  each, desired  output/input  vibration  direction  is  listed.  This 
process  takes  from  1-10  minutes  depending  on  the, number  of  new  choices. 

For  details  see  the  User's  Guide  (Appendix  B). 

The  user  then  activates  the  program  by  typing  "GO  BIODYN",  and  within 
seconds  the  equations  of  motion  are  solved  for  the  desired  transfer  functions 
and  the  results  are  then  printed.  Selected  output  files,  so  Identified  in 
the  CHOICES  File,  can  be  saved  for  later  plotting  or  recall.  From  "GO"  to 
printout,  often  takes  less  time  than  it  did  to' read  this  paragraph.  This  com- 
putational speed  makes  it  feasible  for  a user  to'  sit  at  a terminal  and  itera- 
tively alter  the  INPUT  and  CHOICE  files  to  zero-in  on  the  characteristics  of 
a desired  biodynamic  transfer  function  or  mode. 

The  outputs  can  be  easily  selected  to  provide  any  of  the  several  trans- 
missibility  or  impedance  transfer  functions  called  for  by  the  PIVIB  pilot- 
vibration  performance  program,  mentioned  earlier. 

For  further  manipulation  of  the  resulting  transfer  functions,  such  as 
frequency  response  (Bode)  plots,  calculating  rms  output  in  response  to  an 
input  spectra,  etc.,  the  output  files  can  be  readily  accessed  by  the  user's 
(or  other  Tymshare-resident)  analysis  programs.  A typical  set  is  the  STI 
Unified  Servo  Analysis  Method  (USAM)  programs  which  permit  a variety  of 
systems  computations  and  plots  for  multi loop  pilot -vehicle  systems  of  this 
type  (Ref.  21).  Although  such  extent  programs  are  not  an  intrinsic  part 
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Figure  3.  Procedural  Flow  in  Using  BIODYN-78 
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of  BIODYN- 78 , the  BIODYN  output  files  can  be  easily  transferred  to  them  in 
the  Tymshare  system.  - 

A more  detailed  "User's  Guide  to  BIODYII- 78"  has  been  prepared,  and  it 
is  included  herein  as  Appendix  B (Ref.  2).  This  is  a brief  guide  for  the 
ordinary  user  of  BIODYN,  who  is  somewhat  knowledgeable  in  dynamics,  but  it 
is  not  intended  for  those  interested  in  adapting  BIODYII  to  their  own  computers 
or  in  modifying  it.  For  background  in  operating  the  Tymshare  programs  in 
general,  we  recommend  the  potential  user  read  the  basic  Tymshare  XEXEC  Manual 
(Ref.  22),  as  well  as  the  BIODYN-  78  User's  Guide  (Ref.  2).  Given  a valid 
' • Tymshare/STI  access  agreement,  the  typical  aircraft  engineer  user  can  be 

getting  useful  answers  to  biodyr.amic  interface  problems  within  a few  hours. 

d.  Formal  Demonstration  of  Program 

To  demonstrate  the  simplicity  of  using  BIODYN-78,,  a demonstration  was 
given  on  l6  March  1979  at  AMRL/EM  (WPAFB,  Ohio),  where  a suitable  terminal/ 
printer  and  acoustic  coupler  were  readily  available.  The  demonstration 
was  attended  by  sereral  persons  from  AMRL/BB,  and  a few  from  AFFDL/FGC  (who 
) are  interested  in  using  BIODYN  for  Pilot  Induced  oscillation  problems) . The 

attendees  were  allowed  to  suggest  changes  in  parameters  which  might  be  of 
interest  to  them,  and  these  were  rapidly  computed,  printed  out  and  plotted 
via  a simple  typewriter  plotting  routine  which  is  part  of  STI’s  USAM  program 
(Ref.  21). 

Examples  typical  of  the  demonstrated  cases  are  given  in  Figs.  4 and  5. 
Figure  4 shows-  the  terminal  printout  for  the  standard  test  case,  which  is 
that  used  in  Ref.  17  and  is  explained  in  detail  in  Appendix  B (Ref.  2). 

Here,  the  user  has  selected  a measure  of  head-bobbing j the  head  angular 
acceleration  in  pitch  (8,  rad/s^)  in  response  to  vertical  acceleration  of 
the  platform  (aZp,  m/s2).  The  resulting  transfer  function  factor  files 
are  printed  out  below  and  are  saved  for  later  plotting  (note:  frequencies 
are  in  rad/sec).  Figure  5 shows  the  "quick-look  terminal  Bode  plot"  of 
this  transfer -wherein  the  proximity  of  the  "head-bobbing"  and  "torso  surge" 
modes  at  28-30  rad/s  (=5  Hz)  is  apparent.  This  plot  is  not  currently  part 
of  BIODYN-78  per  se  but  is  available  on  Tymshare  through  STI's  Unified 
. Servo  Analysis  Method  (Ref.  21).  The  time  required  to  obtain  these  print- 

" outs  and  plot  was  about  10  min. 
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List  INPUT  file 


-TYPE  9,19  — 

STIFF  STICK 

1.371  .29  .44  .4124  49. 

1.417  .345  .3  .4152  145. 

.31  2.  13944.  . S57945F.-44  .41  . 61  13944.  1. 

14.  1294.48  29585.27 

2.  1.  84.  44.  2.43123  144.39969  1.  71519.  4. 

5.  .4949  4.  4.  .424437243  .49494949  .48938764 

.8  16.  4.  .S  .45555555  4.  4.  1. 

18.  .15  .8  14.  16.75  544. 

4.14  .3  .439  5.  .44  1.  .4  74.  4. 

13.  .3  5.  .75 

4.  .1  4.  -24.  .25  54.  .45  -34. 

-TYPE  C.Tvp  - ■■  ■ ■ — ■ — List  CHOICES  file 

DTH/DZP 

XXX 

C /AZP 

-CO  BIOOYN  — — '■ Activate  BIODYN 


NOTE:  User  entries  are 
underlined 


enter  INPUT  FILE  NAME:  B.19 

ENTER  CHOICES  FILE  NAME:  C.TNP  } 


Program  prompts 
for  file  names 


ll-Jun-79  21:42 


CASE:  srirr  STICK 


Case  Identifier 


DENOMINATOR 1 


TRANSFER  FUNCTION  COMPONENTS 


• 27577E-17 
( 16.593 

) ( 58.788 

) { 97.341 

) 1 

( 357.45 

) 

((  .26496 

, 14.914 

r 

2.8917  , 

14.524 

((  .24252 

, 12.943 

2.6213 

12.675 

((  .35427 

# 26.761 

9.4846  , 

25.425 

'(  .15112 

, 29.611 

• 

f : 4747 

29.P74 

((  .85762 

, 32.953 

28.261  , 

16.947 

((  .19124 

. 57.483 

• 

14.993  , 

56.422 

( 282.17 


< . 1187  3E*14> 


High  Frequency  gain 
(Denominator) 

First-order  poles 


Second-order  poles 
(£,  in,  (si,  o»y,  _ ,2 

Lev  frequency  gain 
(Denominator) 


NUMERATOR:  DTH/DZP  FILE  NAME?  DTH.TMP 
NEW  FILE 


Program  prompts  for 
file  name 


-. 12319E-1S 
( .44444  ) 

I 95.827  ) 

((  .45587 

((  .97S54E-41, 
((  .84471  , 

II  .23595  , 

((-.26361 
<-.29394E»«7> 


( .44844 

) ( 16.594 

) ( 

18.284  ) | . 

( 232.81 

) ( 498.67 

) 

1 

14.396 

, 4.7392  , 

9.2528 

>>  > 

16.243 

, 1.5846  , 

16.125 

)) 

33.835 

, 2.  . 445.  , 

18.321 

)>  

51.815 

, 12.226  , 

54.352 

284.74 

,-75.464 

274.66 

))  J 

High  frequency  gain 
(Numeratcr)' 
First-order  zeros 


( Second-order  zeros) 


Liv  frequency  gain 
(ilumerat-or)' 


NUMERATOR:  C /AZP 

. 38429E-14 


( 11.444  ) 

( 14.817 

) { 22.254 

) ( 

23.472  ) | 

(-38.117  ) 

( 184.44 

) 

V S 

imi’ar  to 

((  .19965  , 

12.356 

, 2.4668 

12.147 

>>  r — s 

Dove 

((  .85427  , 

17.818 

, 15.222  , 

9.2626 

>>  1 

((  .14663  . , 

28.944 

, 4.2381  , 

28.591 

)»  1 

((  .97441  , 

147.45 

, 143.84  , 

26.424 

>>  ) 

<-.57327E»14> 

TRU  « 17.42 

ll-Jun-79 

. 21:43 

EXIT 


Figure  k.  Typical  BTPDYN-78  Computer  Printout 
for  Example  Cast 
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-GO  (Y7LI8)  USAM2 


ll-JUN-79  21:47 


WHAT  NOW  ? FIN 
FILE  NAME  ? DTH.TMP 


User  entries  are  underlined: 

Selects  file-name 
for  plotting 


WHAT  NOW  ? BPL 
ENTER  WO  WF:  1. , 133. 


Requests  "Bode-Plot'* 
end  sets  frequency  range 


AMPLITUDE, + - DB 
23.  9.  -29. 


-8a.  -laa. 


iaa.aaa> 


79.829> 


53. 103> 
39.899> 


28.2 aa>  + 


• 29.399> 


— i4.iaa> 
4 

M n.aaa> 


•*•  7.aaa> 


s.aia> 

3.983> 


2.8ia> 


2 . aaa> 


i.4ia> 


i.aaa> 


+ * 

+ * 

+ • 

* 

*+ 

* ♦ 

* * 

* + 

* ♦ 


■S  5 

'i 

09  U 

4>  o) 

l S 

> * 

(0  © 

9 5 


+»  . s 
22  ? 
M h e 

2®  a 

2 f.  * 

1°  5 

& X tc 
O O -H 

£ S'  S 


233.  iaa. 

PHASE,*  - DEG 


-203.  -333. 


Figure  5.  Typical  USAM2  Terminal  Bode  Plot  for 
Check  Case 
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standard  I&E  SEKI-I/XiARITHKIC  3 cycles  x 70 
divisions  paper. 


The  formal  demonstration  and  the  attached  User's  Guide  (Appendix  B) 
completes  the  primary  development  of  our  biomechanical  model  for  feed- 
through of  vibrations  to  a pilot's  hands  and  head. 

C.  CONCLUSIONS 

Through  the  use  of  modern  computer  techniques  and  closed-loop  neuro- 
muscular models.,  the  biomechanical  feedthrough  of  seat  vibration  of  the 
various  parts  of  the  torso,  head,  arm  and  hands  (on  a restrained  ce-trol 
stick)  has  been  fairly  accurately  and  efficiently  modeled. 

The  situation  requires  a model  having  80-odd  postural  and  physiological 
parameters  to  define  properly  the  general  case  of  vertical  and/or  fore-aft 
vibrations.  An  internal  routine  automatically  trims  the  Tor so/Head /Limb 
system  at  the  desired  ingles  and  computes  linearized  equations,  which 
are  solved  for  their  eigenvalues  and  resulting  transfer  functions  in 
seconds.  In  distinction  from  most  past  models,  this  model  represents 
simultaneously  several  degrees-of- freedom  e.g.,  vertical.  shouJder  motion, 
angular  head-bobbing  notion,  constrained  elbow  motion,  stick-grip-interface 
force,  and  stick  angular  motion,  among  many  others.  These  features  make 
the  BIODYN-78  model  well  suited  for:  routine  computations  of  biomechanical 
transmissibilities,  (e.g.  those  required  as  inputs  to  the  PIVTB  Program  of 
Ref.  20),  matching  of  test  data,  systematic  analysis  of  variations  in  body 
parameters  (e.g.  to  explore  population  extremes),  optimization  of  control  , 
stick  "feel"  parameters  (i.e.,  to  understand  or  avoid  Pilot  Induced  Oscil- 
lation), or  to  compute  control  transients  due  to  accelerations  imposed  on 
the  pilot  by  external  forces. 

Table  2 compares  the  final  results  reported  here,  as  well  as  in  .-'ior 
Interim  Scientific  Reports,  with  the  criteria  imposed  on  page  i. 

This  fifth-year  effort  was  focused  on  the  last  objective.  The  User's 
Guide  (included  here  as  Appendix  B),  and  the  formal  demonstration  of 
BIODYN-78  operation  using  a remote-access  terminal  in  the  Aerospace  Medical 
Research  Laboratory  (at  WPAFB,  Ohio),  completely  meet  that  objective.  The 
other  objective  of  the  FY  1978  program  was  to  evolve  a model  for  the  effects 
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TABLE  2.  ASSESSMENT  AGAINST  AFOoR  CRITERIA  PJR  ENGINEERING  MODELS 


CRITERION 

a.  Based  on  sound  empirical 
data 


b.  Capture  essential  phenomena  of* 
interest  to  users,  without 
excessive  restrictions  or 
details . 


c.  Compatible  with  more  detailed 
physiological  models.' 

i 


d.  Put  in  fora  usable  by  others. 


> 

3R-1037-I 


ASSESSMENT 

• Model  has  been  fit  to  several 
comprehensive  sets  of  quasi- 
random vibration  data. 

• Can  cope  with  all  likely  pilot 
postures,  including  supine. 

• Can  handle  problems  of  various 
kinds,  of  stick,  arm  rests,  torso 
restraints  (linearized). 

• Shows  common  resonant  body  modes 
seen  in  several'  sets  of  wideband 
data  (0.5  to  20  Hz) . 

• Shows  dominant  modes  at  limbs, 
shoulder  and  head/neck,  but  lacks 
detail  for  rib-cage,  visceral  and 
pelvic  motions. 

• Needs  comparison  against  more 
recent  detailed  models  of  Belytschko, 
Kaleps,  etc. 

• Is  now  in  a national  time-shared 
computer  service,  accessible  by 
standard  remote  telephone  terminals. 

• Users  Guide  and  test  case  available 
to  any  potential  user. 

• Access  procedures  need  to  be  made 
more  readily  available  to  general 
users . 
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of  head  vertical  plane  vibrations  on  eye  and. image  motions;  (as  influenced 
by. the  vestibulo-ocular  and  fixation  reflexes).  This  was  met  by  the  model 
detailed  in  the  previous  I3S  (Ref.  19)  and  summarised  in  Appendix  A here. 


In  summary,  nearly  all  of  the  goals  of  this  5-year  program  have  been 
met,  with  some  re-direction  away  from  detailed  lateral  vibration  models 
(ours  remains  the  simple  one  we  presented  ir.  Ref.  25)  towards  the  more 
detailed  head  and  eye  motions,  noted  above.  This  latter  work  has  not  yet 
been  incorporated  into  the  computerized  BIODYN-73  program.  Since  it  seems 
to  explain  seme  of  the  Image  motion  and  visual  acuity  differences  between 
vibrating  targets  and  vibrating  heads,,  (see  Appendix  A)  this  link  promises 
to  complete  the  chain  of  dynamic  elements  between  seat  vibration  inputs 
and  visual  performance  decrements,  including  head-mounted-displays. 

D.  RECOMMENDATION 

Even  though  the  BIODYN -78  model  is  now  in  a readily  usable,  computer- 
ized form;  complete  with  a User's  Guide;  some  further  improvements  are 
needed: 

a.  The  ease  of  access  of  general  users  to  the  Tymshare  resident 
BIODYN-78  program  must  be  improved.  At  present,  a user  must 
contact  STI  for  a special  users  agreement  (at  no  profit  to 
STI),  ir.  order  to  gain  access  to  the  program.  Having  STI 
assist  new  users  in  applying  BIODYN-78  is  fine,  but  the  pro- 
cedure is  awkward  for  military  facilities  and.  some  large 
companies  with  their  own  in-house  computer  facilities.  We 
recommend  that  some  means  be,  worked  out  for  Tymshare  support 
of  BIODYN,,  possibly  with  an  optional  STI  advisory  role. 

b.  The;  present  program  only  prints  out  or  files  the  requested 
transfer  functions  in  pole-zero  numerical  format.  The 
addition  of  a quick-look  computer  Bode-plotting  routine. 

(such  as  demonstrated  herein  via  our  USAM  program),  would 
greatly  facilitate  the  use  of  BIODYN  in  routine  experimental 
and  engineering  development  projects. 

c.  The  head-eye  model  should  be  incorporated  into  BIODYN  for  use 
in. analyzing  aircrew  visual  performance  under  vibration. 

d.  Efforts  to  publicize  the  availability  of  BIODYN-78  to 
interested  researchers  and  development  engineers  should 
be  undertaken,  to  assure  that  the  whole  program  reaches 
its  potential  usefulness  to  the  aerospace  and  clinical 
community. 
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progress  nr  measurer;  and  jcbelisg  the  effects  or 

LOW  FREQUENCE  V 13 RATION  ON  PERFORMANCE 
by 

Henry  R.  Jex  and  Raymond  E.  Magdaleno 
Systems  Technology,  Inc. 

1 57f6  S.  Hawthorne  Boulevard 
Hawthorne,  California  902fO 
U.S.A. 


StDMARY 

Several  facets  of  the  comprehensive  blodynaaic  modeling  program  presented  at  the  AGAPD  Aerospace 
Medical  Panel  Meeting  at  Oslo,  '97^,  have  been  successfully  completed  and  are  reported  here.  The  objec- 
tives and  approach  (an  ensemble  of  physical  models  of  Just-adequate  complexity,  with  adjustable  or  adap- 
tive parameters)  are  first  reviewed. 

The  development  of  a variety  of  lumped  parameter  models  to  explain  and  codify  the  known  data  on  low- 
frequency  vibration  effects  and  to  predict  likely  effects  in  new  situations  has  been  brought  to  a useful 
level.  These' are  described  with  presentations  of  typical  validation  data: 

a)  BIDDYN-T8,  a user-interactive  computer  program  which  is  capable  of  modeling  the  biomechani- 
cal properties  of  a variety  of  pilot/crevaan  posture  and  control  situations  in  which  active 
neuromuscular  systems  are  involved. 

b)  PIVIB,  a user-instructed  batch  program  (developed  by  others)  for  modeling  effects  of  vibra- 
tion on  tracking  performance,  and  which  requires  information  supplied  by  BIODl'N-78  in  its 
us*. ' 

The  relationship  and  applications  of  these  and  other  related  new  models  are  discussed  with  respect  to 
their  development  status  and  pc tential  applications. 

Specific  recommendations  are  made  for  more  refined  experimental  data  (e.g.,  simultaneous  accelera- 
tions on  various  body  locations  and  better  postural  and  dynamic  mode  shapes  via  cinematography,  etc.)  and 
interface  compatibility  among  various  models. 

OBJECTIVES 

This  is  a progress  report  on  recent  work  and  achievements  in  measuring  and  modeling  the  effects  of 
low- frequency  vibration  (0.2-20  Hz)  on  the  performance  (visual,  tracking,  etc.)  of  vehicular  crewmm. 
Although  the  main  research'  relates  to  pilots  and  operators  of  equipment  In  aircraft,  there  are  examples 
from  high-speed  ships,  and  the  resulting  technology  certainly  applies  to  land  vehicles. 

Typical  problems  for  which  these  models  are  being  developed  are  shown  in  Fig.  1.  Vibration  effects 
on  performance  usually  occur  at  the  operator  interfaces  with  his  environment,  such  as:  the  seating/ 
restraint  system  (and  any  isolation  features  such  as  cushions  or  springs);  displays  (whether  or  not  they 


EYE' HEAP/ *IA«ET  DYNAMICS: . 

• Analyze  vibration  nterference  with  vision 

• Biodynamic  design  ol  helmet-mounted  display! 
under  high  Gr  ard  buffet 


RESTRAINT' SEATING!  AND  ISOLATION  SYSTEMS: 

« , Anal  yze  performance  consequences  d competing 
restraints  end  cushions 

s Design  new  restrai  nt  and  isolation  systems  ■ 

s Analyse  perforrrarce  effects  d various 
postural  situations  U.*..  supine  seating) 

e Ride  quality  impiic  rtwns  on  various  tests 
end  locations  met rcrett  . 


DISPUYS: 

a Ana  tyre  display  mounting,  isolation, 
stabilisation  concepts 

a Analyse  effects  ot  size,  viewing  distance, 
orientation,  etc 

CONTROIS'MOTOR  COORDINATION: 

, a Analyze  competing  control  configurations 
la.  g. , fly-by-wire  vs.  canter  stick!  lor 
la)  affects  of  rough  air  eaeited  bending 
modes  for  large  aircraft,  or  Ibl  high  g 
buffet  for  fighters 

a Design  vibration  resistant  controls 


a Optimize  among  various  tradeoffs  prior  to 
testing,  also  simulations  to  operational  ■ 
situations 
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figure  1.  Problem  Area*  of  Interest 
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move  with  the  seat  or  cockpit);  controls  (especially  tracking  control  sticks  or  devices  held  and  aimed  by 
the  operator)!  eye/head  motions  as  they  affect  the  perceived  image  motions  (and  effect  of  helmets  and 
helmet-mounted  sights).  Finally,  there  are  the  overall  tradeoffs  between  better  seating  and  restraints 
versus  better  controls  vers--,  better  displays  (e.g. , vibration-compensated,  displays,  head-mounted  dis- 
plays, etc.).  For  such  problems  as  these,  one  needs  "uniformly  valid  first  approximation”  models,  which 
can  capture  the  essential  features  of  a variety  of  vibration  Interface  problems  in  a computationally  effi- 
cient manner.  That  is  tne  overall  goal  of  this  research. 


APPROACH 

Overall 

A more  detailed  discussion  of  these  problems  and  possible  approaches  toward  solving  them  was  pre- 
sented at  the,  Aerospace  .'Medical  Panel  Specialist's  on  "Vibration  and  Combined  Stresses  in  Advanced  Sys- 
tems," at  Cslo,  Norway,  in  April  1971*  (°ef.  1).  After  reviewing  the  relevant  definitions  and  measures  of 
task  performance,  behavior,  stresses,  and  strains,  that  paper  considered  three  approaches: 

• A purely  experimental  Investigation  of  all  possible  situations  (which  would  never  be  even 
partially  completed). 

• An  ad  hoc  simulation  of  each  problem  as  it  ar'  se  (which  gives  relevant  data  for  the  case 
at  hand  but  cannot  be  extrapolated  to  new  situations). 

• A comprehensive,  empirically  based  theory  (which  requires  carefully  coordinated  cycles  of 
modeling,  validation/refinement  experiments*  and  the  development  of  a user-oriented  "cata- 
log" of  models  and  data  that  can  be  extrapolated  to  new  situations  even  before  a proper 

ad  hoc  simulation  can  be  developed). 

The  comprehensive  theory  approach  waa  recommended  in  Ref.  1 for  a number  of  reasons,  given  therein, 
and  the  USAF  Office  of  Scientific  Research  and  Aerospace  Medical  Research  Laboratory  have  sponsored  the 
majority  of  the  program.  This  paper  reports  progress  as  of  mid  1978,  at  which  time  the  first  really 
usable,  albeit  embryonic  models  are  available. 

Scops  and  Status 

As  noted  above  and  discussed  In  Refs.  1-3,  the  basic  approach  is  to  use  simplified  physical  models 
with  dynamic  elements  analogous  to  those  of  the  body/liab/heaf/eye  system  and  their  interfaces  with  the 
sept,  display,  and  controls.  The  complete  three-dimensional  assesblage  of  bones,  muscles,  and  organs, 
plus  the  complex  neuromuscular  control  system,  is  far  too  complex  to  model  efficiently  enough  for  prac- 
tical solutions  of  the  problems  described,  and  such  complexity  is  not  really  needed.  Instead,  only  those 
elements  required  for  a likely  class  of  problems,  such  as  those  of  Fig.  1,  are  employed. 

Most  of  this  work  to  date,  has  emphasized  vertical  and  fore/aft  vibration  (e.g. , Refs,  b-6),  with 
some  work  on  lateral  vibration  (e.g..  Refs.  2 and  7)  and  very  Halted  work  on  combined  vibrations  (e.g.. 
Ref.  8). 

Figure  2 shows  the  fundamental  model  structure  used  for  evaluating  biodynaraic  interference  with 
visual  and  control  tasks,  adapted  from  Ref.  1.  The  main  blocks  delineate  the  "p>man  operator,"  "con- 
trolled element"  (which  is  also  the  source  of  vibrations  vis  the  disturbance  Inputs),  and  "manipulated 
object^  (such  as  a book,  aimed  device,  or  simply  a visually  fixated  external  object.  The  main  interfaces 
at  the  display,  seat,  and  controls  axe  also  shown.  Within  each  major  block  are  the  main  elements  ve  have 
modeled,  e.g.,  seat-interface,  body,  head,  eyes,  limb,  etc.  Ws  will  expand  on  these  later. 

The  checks  (,/)  or  question  marks  (?)  beside  each  element  denote  the  status  of  its  math  model.  The 
earlier  ( 197'1 ) status  Is  shown  dashed  to  reveal  where  progress  has  been  made.  There  are  few  left, 
implying  that,  the  models  «xe  approaching  a new  level  of  maturity  in  which  they  are  ready  for  more  wide- 
spread validation,  refinement,  and  application. 

The  selection  of  specific  models  from  the  available  ensemble  starts  with  Fig.  2,  using  Its  blocks  to 
define  the  task  and  Its  measures  of  performance!  the  controlled  element  with  its  dynamic  elements  includ- 
ing structural  modes,  seat-responss  properties,  and  control  excitation  of  vehicle  motions;  the  forcing 
functions  (commands  to  bs  followed  and  disturbances  to  be  suppressed  or  which  produce  vibration);  the 
interface  properties  such  as  display  distances  and  sizes;  control  stick  "feel"  properties;  etc.  The 
specification  of  all  the  dynamic  elements  of  Fig.  2 is  essential  to  any  application  and  requires  a fairly 
complex  model  structure  and  an  input-descriptor  file  containing  several  dozen  parameters. 

The  next  step  is  the  calculation  of  the  biodynaalc  "feedthrough"  or  transmissibility  to  the  seat, 
head,  eyes,  and  image  motions  (for  visual  performance  effects)  and/or  to  limb  and  control  stick  (for 
minuel  control  performance).  All  of  the  models  used  here  have  special,  features  to  allow  the  redundant  , 
and  active  neuromuscular  forces  due  po  body/limb  coupling  with  the  control  sticks  or  arm  rests  to  be 
efficiently  modeled.  The  manual  control  situation  precludes  the  use  of  the  pasalve  "open  chain"  anthropo- 
morphic models  now  highly  developed  for  crash  research  (e.g.,  Ref.  9>. 

The, model  and  user's  program  ve  have  been  developing  to  coepute  these  biomechanical  transmissiblli- 
ties  is  called  "3l09Ytl-78"  (BIODYHaoics,  and.  the  number  denotes  the  model  year),  and  its  scope  and  role 
are  depicted  on  Fig.  3,  which  la  an  expansion  of  the  center  region  of  Fig.  2 (Ref.  10).  Even  for  rela- 
tively simple  situations  involving  tracking  a CRT  display  with  « spring- restrained  control  stick,  thers 
axe  numerous  biodynaalc  elements  In  the  causal  chain,  and  dozen*  mace  parameters  must  be  specified  to 
define  the  operator's  size,  posture,  tenseness  (neuromuscular  gains),  eye-to-display  locations,  etc.  A* 
noted  by  the  numbers  in  parentheses  in  Fig . 3,  there  are  • seat/pelvis/interface  parameters  ( such  as 
seat  fores  gradient,  damping,  lower  body  effective  mass,  compliance  and  damping) , 10  torso  parameters. 
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21  neuromuscular  (TIM)  items,  1?  limb./stiok  interface  or  arm  rest  parameters,  and  il*  head/neck/shoulder 
parameters.  These  are  in  the  1978  version,  which  computes  limb/stick  transaissibility  or  forte  impe- 
dance as  well  as  head-point-of-regard  and  translation,  shoulder  motions,  etc.,  as  illustrated  in  Fig.  l>. 
Only  the  main  body  links  are  shown  in  Fig.  la.  > The  basic  equations  of  motion  are  nonlinear  in  the  pos- 
tural and  neuromuscular  w-riables,  but  the  equations  are  automatically  trimmed  about  the  average  operat- 
ing point  and  linearized  t.  'ations  are  used  to  compute  the  vibration  effects. 


A first  cut  at  a model  for  vibration- induced'  image  motions  has  been  made.  (Kef.  11)  which  takes  as 
inputs  the  head  rotations  and  translations  and  accounts  for  the  vestibulo-ocular  reflex  (which  roughly 
stabilizes  the  eye  Inert tally  up  to  about  6-3  Hz)  and  the  fixation  reflex  (which,  can  follow  display  image 
notions  up  to  about  1-2  Hz).  It  is  denoted  on  Fig.  J as  "BIODYN-79"  in  anticipation  of  its  inclusion  in 
the  current  3X0DY3'79  computer  program,  and  is  illustrated  in  Fig.  ?.  About  It  eye/head  quantities  are 
specified  along  with. five  or  six  display  interface  parameters.  (For  details  see  Ref.  11-)  We  also  plan 
to  include  the  additional  effects  of  a helmet  and  helmet -mounted  displays  or  starts  in  the  1979  version. 

As  the  third  step  the  computed  vibration-induced  motions  of  the  head,  stick,  3houlder,  and  points  of 
regard  (at  the  display)  of  the  head  and  eye  are  used  to  determine  the  visual  task,  performance  decrements 
via  empir'cai.  relationships  between  these  motions  and  perceptual  or  motor  performance  (e.g.,  detection 
errors,  tracking  "indifference  threshold,"  etc.).  The  connection  between  the  net' performance  of  various 
tasks  and  the  computed  (or  actual  limb  or  head)  motions  remains  the  weakest  link  in  the  chain  from  vibra- 
tion input  to  measured  task  performance. 

Further  systematic  experiments  of  the  type  originally  performed  by  Drazin  (Ref.  12)  and  lately  by 
Benson  and  Barnes  (Refs.  13  and  ib)  need  to  be  performed  to  provide  a data  base  for  this  final  link  in 
the  visual  performance  chain. 
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Figure  5.  Model  for  Image  Motions  Due  to  Head  Vibration 

When  tracking  performance  (or  performance  on  a manipulated  object  task)  la  desired,  the  fourth,  and 
most  difficult,  step  is  necessary.  All  of  the  vehicle  dynamics  and  biomechanical  transmittar.ces  must  be 
included  In  a complex  model  of  the  human  operator's  adaptive  tracking  strategies.  The  "classical"  models 
for  doing  this  are  typified  by  those  of  McRuer  and  his  associates  'e.g.,  Refs.  If  and  16).  They  use 
lumped-parameter  submodels  for  the  perceptual,  equalization,  and  r.  ui  ^muscular  processes,  with  adjust- 
able parameters  which  covary  with  various  task  variables,  including  vibration  (e.g.,  see  Ref.  2).  A 
current  version  of  this,  the  so-called  "Compleat"  Human  Tracker  model  is  summarized  in  Fig.  6. 

The  modern  equivalent  to  this  is  the  optimal  control  theory  model  of  Kleinman,  et  al;  (e.g.,  Ref.  17), 
in  which  the  catalog  and  adaptation  "rules"  of  classical  human  operator  model  are  replaced  by  an 
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optimal  control  law  algorithm  whose  cost  functional  is  adjusted  (based  c.n  fitting  empirical  data)  to  pro- 
vide the  covaryir.g  parameter  trends  observed  in  experiment::. 

The  optimal  control  model  for  tracking  has  been  highly  refined  and  applied  to  modeling  performance 
under  various  axes  and  types  of  '-lbration  by  Levison  (e.g.,  Refs.  6-3)  and  has  been  put  in  a user-oriented 
computer  program  as  "FIVIB"  (Pilot  VIBration,  Re.f.  18),  which  in  accordance  with  our  dating  scheme  will 
be  designate!  herein  as  "PrVIB-77."  As  implied  on  Fig. . '3,  PIVIB-77  can  use  BIODYN-78  (or  experimental 
data,  if  available)  to  provide  the  various  biomechanical  transmissibilities  required.  The  current  PIVI3 
has  some  oversimplified  assumptions  regarding  the  performance  decrements  due  to  head-point-of-regard  motion 
(e.g.,  It  does  not  yet  properly  account  for  vestibular  or  image  tracking  effects  of  the  eye's  control  sys- 
tem). Nevertheless,  it  does  give  a good  first  approximation  to  performance  changes  for  a wide  range  of 
task  variables  (such  as  input  spectra,  stick  scaling,  vibration-induced  blur  threshold).  As  should  be 
expected  from  the  discussion  of  Figs.  2 and  5,  hundreds  of  parameters  are  required  in  the  PIVIB  input 
files,  because  it  has  to  cover  all  of  the  elements  in  Fig.  2. 

We  have  shown  that  each  of  the  four  steps  in  assessing  vibration  effects  on  performance  now  has  math 
models  of  varying  degrees  of  confidence  and  ease  of  use.  The  selection  of  individual  models  depends  on 
the  axes  of  vibration  and  of  manual  control,  the  type  of  vibration,  and  the  availability  of  purely  empiri- 
cal data  (e.g.,  experimental  frequency  responses)  or  the  need  for  prediction  before  the  experiment  or 
design  is  frozen. 


VALIDATIONS 

The  biomechanical  model,  BIODYlt-78,  and  its  antecedents  have  been  progressively  refined  and  validated 
as  experimental  data  of  suitable  format  and  quality  became  available.  Typical,  fits  of  the  body,  limb,  and 
head-point-of-regard  data  are  shown  in  Fig.  Jc  from  one  of  a series  of  experiments  with  quasi-random 
vibration  (i.e.,  the  sum  of  five  sinusoids  of  non-simple-harmonic  ratios).  Further  details  are  given  in 
Ref.  3-  The  main  point  to  note  from  Fig.  kc  is  that  one  model  (albeit  complex)  simultaneously  fits  all 
of  the  several  degrees  of  freedom  measured.  More  experiments  are  needed  with  simultaneous  measurements 
at  several  well-defined  body  points,  e.g.,  shoulder,  sternum,  head  (2  places),  and  arm  bones  or  joints, 
along  with  the  main  seat  input  (e.g.,  via  a 10  cm  "ISO”  disk  with  accelerometers  between  the  outtocks  and 
seat) , Since  small  differences  between  large  accelerations  are  important,  very  high  grade  instrumenta- 
tion and  inputs  are  required. 

Figure  -b  shows  another  example,  also  described  in  Ref.  3,  where  wideband  vehicle  vertical  motions 
In  the  0.;  to  5 Hz  range  (from  the  high-speed  ship  simulation  program  described  at  Oslo  in  Ref.  1)  were 
imposed  on  crewmen  sitting  in  various  degrees  of  postural  relaxation  (here  relaxed,  eyes  closed).  The 
main  point  is  that  the  multimodal  character  of  the  data  is  well  matched  by  the  model.  One  subject  noted 
that  he  "tilted  his  hips  forward"  (dashed  line  in  Fig.  kb)  to  relieve  some  of  the  head  bobbing;  the  model 
correctly  captures  this  effect. 

The  head/eye/control  image  motion  model  of  Fig.  5 is  much  less  well  validated,  but  it  shows  promise. 
It  relies  on. the  fact  that  under  moat  low-frequency  vibrations  the  eye  only  moves  a degree  or  two  (not  ' 
enough  to  trigger  visual  saccades)  and  thereby  permits  a quasi-linear  modeling. 

At  the  bottom  of  Fig.  7 are  some  image  motion  prediction*  made  for  three  types  of  display /head  inter- 
face; target  alone  moved  (head  fixed);  target  fixed  (head  and  eye  moved);  and  a head-mounted  display 
(image  moving  with  the  head).  Notice  the  different  bandvldths  of  image  (error)  motion,  as  mentioned  ear- 
lier — poor  for  target  moved  and  good  for  head  moved.  More  interestingly,  the  head-mounted  display  gives 
rise  to  large  image  motion  errors  because  the  counterpitching  tendency  of  the  eye*  works  against  a head- 
mounted  display.  Vibration-induced  decrements  in  digit  reeding  accuracy  a*. measured  recently  by  Benson 
and  Barnes'  experiments,  are  shown  at  the  top  of  Fig.  7-  The  similarity  of  trends  in  computed,  image 
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Figure  7.  Computed  Image  Motions  and  Digit  Reading  Performance 


motion  and  measured  accuracy  suggests  some  validation  for  assuming  that  visual  performance  can  be  related 
to  image  motion.  Again,  muph  more  experimentation  is  needed,  and  we  recommend  that  eye  motion  be  measured 
in  as  many  cases  as  feasible. 


The  overall  ability  of  the  complex  PIVIB-77  model  structure  and  algorithms  to  model  manual  control 
performance  is  demonstrated  in  Fig.  8,  which  includes  overall  error  scores  (top)  and  control  actions 
(bottom)  from  three  different  experiments  performed  at  the  USAF  Aerospace  Medical  Research  Laboratory. 

A couple  of  generil  points  are  illustrated  by  Fig.  3: 

• lhe  PIVTB  model  (after  adjustment  of  its  algorithms)  was  able  to  fit  all  three  experiments 
fairly  well  without  further  changes,  both  for  error  and  control  action. 

• The  effects  of  vibration  versus  static  on  task  error  are  generally  nsc  very  large  and  are 
often  within  the  range  of  individual  subjects.  Adversely  sensitive  control  stick  gains  or 
postures  can  greatly  increase  the  vibration  effects,  espe.'-tally  on  the  rms  control  activity, 
as  shown. 

• Big  effects  occur  due  to  input  scaling,  stick  and  display  sealing,  etc.  These  effects  are 
. currently  modeled  by  PXVI3.  ’ 


, PTVIB  is  hot  computationally  very  efficient,  because  it  incorporates  optimum  Kalman  filter  subrou- 
tines which  require  tbs  solution  of  Riccati  equations.  However,  it  has  been  put  in  usable  form  with  a 
user's  guide,  ar-i  it  is  ready  to  be  tested  and  further  refined. 
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An  objective  of  these  models  is  to  permit  prediction  of  problem  areas  or  solutions  before  actual 
simulation  data  are  available.  One  such  prediction  of  BIODYH-78  is  shewn  in  Fig.  9,  which  compares  the 
he*d-line-of-sight  traiismissibillty  (which  the  eyes  have  to  compensate  for)  for  a standard  15  deg  tilt- 
back  fighter  aircraft  seat  with  an  advanced  6?  deg  tilt-back  seat  designed  for  High  Acceleration  Cockpits 
(HAC)  (for  details  see  Ref.  J).  An  order  of  magnitude  increased  sensitivity  to  vibration -induced  head 
bobbing  is  predicted  if  the  head  or  neck  is  not  supported  (as  assumed  herel.  Vibration  experiments  on  • 
such  configurations  seem  warranted,  with  comprehensive  measuremente  of  seat  and  torso  motion,  head  and 
limb  motions,  and  tracking  and  visual  identification  task  performance. 

Other  predictions  have  been  made  as  to  the  influence  of  an  arm  rest  oh  vibration  feedthrough  (depends 
on  the  percentage  of  arm  weight  carried  by  the  rest  and  on  the  limb  angles  Involved)  and  on  the  effects  of 
various  postures. 

The  usefulness  of  such  models  for  this  sort  of  prediction,  parameter-tradeoff  study  depends  on  how 
convenient  is  its  access  and  how  easy  it  is- to  modify  the  parameters.  BIODYN-"^  has  been  designed  for 
user-interactive  operations  on  the .Tymshare,  Inc;,  tlmeshared  computer  network,  and  it  could  eventually 
be  accessed  by  anyone  with  a simple  teletype  terminal  with .telephone  coupler  to  Tymshare  (not  yet  avail- 
able,' however).  Among  the  advantages  are  the  very  convenient,  well  documented,  and  powerful  Editor  arid 
Executive  modes  for  Tymshare,  vhich  make  it  easy  to  change  the  variables,  select  new  inputs  or  outputs. 
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Figure  9.  Use  of  Model  for  Problem  Prediction, 
(increased  Head  Motions  in  High-G  Seats) 


and  store  and  plot  the  results  on  a simple  terminal.  We  hope  that  this  scheme  can  be  made  more  available 
for  general  use  and  we  encourage  others  to  do  the  same.  PIVIB  is  a batch-oriented  program  presently 
available  on  the  WPAFB  CDC-6600  computer  for  USAF  projects  only. 


PROBLEMS 

■ ‘(either  BICDYN-78  nor  PIVIB-T7  are  complete  programs,  and  numerous  gaps  and  weak  areas  exist  which 
need  attention.  Sotos  of  .these  are  as  follows: 

1)  BIODYN-78  needs  to  be  extended  to  incorporate  the  lateral  model  previously  derived  in 
Ref.  2,  the  head-eye  model  of  Ref.  11,  and  helmets  and  helmet-mounted' displays  and  active 
neuromuscular  control  of  the  head/neck  system. 

2)  PIVIB-77  needs  better  interface  models  with  respect  to  image  motion  effects,  free  stick 
effects,  as  well  as  a more  convenient  interactive  mode  of  operation  (in  addition  to  the 
efficient  batch  mode). 

3)  The  links  between  head,  hand,  and  eye  motions  and  visual  perception  effects  need  to  be 
forged  on  a sound  scientific  basis.  The  theory  should  be  used  to  guide  the  experimental 
plans  and  measurements. 

*• ) The  validity  of  the  neck-head-vestibular-eye  system  needs  to  be  extended  to  very  low  fre- 
quencies (below  0-5  Hz)  where  some  anomalies  have  consistently  shown  up  (e  g , see  Fig.  Lb). 
Eye  movement  data  are  essential  here,  and  postural  controls  must  be  imposed. 

’ ) More  users  have  to  try  to  apply  these  models  to  their  data.  One  ideal  candidate  is  Row- 
lands' extensive  compilation  of  seat  to  head  and  shoulder  transmissiblLitles  under  a wide 
variety  of  postures  (Ref.  19).  The  main  obstacle  here  is  the  lack  of  documented  postural 
information  (e.g.,  two-view  photos)  or  of  vibrational  "mode  shapes"  at  the  salient  peaks 
and  dips  of  each  frequency  response,  (e  g.,  via  high-speed  cinematography).  Please  docu- 
ment your  postures  and  dimensions  for  each  subject; 

' 6)  The  movement  of  the  eyes  under  vibration  is  very  difficult  to  measure,  because  nothing  of 
any  mass  can  be  attached  to  the  face  ( e lectro-oculograas  are  no  good  when  facial  tissues  ' 
Jiggle  around)  and  the  motions  are  very  small.  The  TV  scanning  oculometer  such  as  Honey- 
well's would  seem  ideal,  but  its  effective  sampling  lag  Is  excessive  for  vibration  work.  . 
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"USER'S  GUIDE  TO  BIODYN-78,  AN  INTERACTIVE  COMPUTER 
PROGRAM  FOR  MODELING  BIODYNAMIC  FEEDTHROUGH 
TO  A PILOT'S  HANDS  AND  HEAD" 


A.  INTRODUCTION 

BIODYN-78  is  a self-contained  program  for  calculating  transmissibili- 
ties  (transfer  functions)  between  vibration  inputs  (vertical  and/or 
fore-aft)  and  various  biodynamic  outputs  for  a seated  pilot  gripping  a 
center-  or  side-stick  and  viewing  a display  while  (possibly)  engaged  in 
a tracking  task.  Model  rationale, ' development,  and  validation  were 
reported  in  Ref.  1 , included  here  as  an  appendix. 

The  purpose  of  this  guide  is  to  briefly  review  and  update  the  model's 
geometric  and  functional  parameters,  and  to  describe  the  two  files  needed 
to  run  BIODYN-78  (model  parameter  "INPUT"  and  transfer  function  "CHOICES"). 
A "test  case"  example  execution,  using  the  parameters  from  one  of  the  vali- 
dating experiments  (Ref.  l),  is  used  to  illustrate  typical  outputs. 

The  BIODYN-78  program  was  developed  on  the  Tymshare,  Inc.  timesharing 
* computer  network  (USA  and  Europe) . The  reasons  for  having  BIC.DYN-78  resi- 

dent on  such  a facility  are: 

• Ease  of  access  by  a simple  typewriter  terminal  with 
modem,  without  "batch"  delays.  (The  program  runs 
each  case  in  seconds.) 

• Ease  of  creating,  modifying,  storing,  and  manipula- 
ting files  of  input  and  output  data  using  the  resi- 
dent EDITOR  program. 

• Availability  of  excellent  Users'  Manuals  to  cover 
the  basic  procedures  in  using  a given  facility. 

This  is  especially  important  because  BIODYN  is 
intended  for  routine  use  by  aircraft  system  design 
engineers  who  often  are  not  computer  ekperts. 

In  view  of  the  last  point  above,  we  will  provide  minimal  information  on  the 
basic  computer  procedures  and  file  operations,  and  will  assume  that  the 
user  is  familiar  with  Tymshare's  XEXEC  manual  (Ref.  2). 

I 
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While  BIOBYN-78  resides  in  Tymshare,  Inc.'s  POP- 10  based  timesharing 
network,  it  has  also  been  recently  test  run  on  Control  Data  Corp.  's 
CYBERNET  NOS- 175  Service  (interactive  mode),  with  identical  ease  and 
results.  If  there  are  any  questions  on  the  use  of  BIODYN-73  on  these 
facilities,  please  contact  either  of  the  authors  at: 

Systems  Technology,  Inc. 

13766  S.  Hawthorne  Blvd. 

Hawthorne,  CA  90250 

(21 3)  679-2281 

A simplified  flow  diagram  for  the  procedures  used  in  applying  BIODYN-78 
is  given  Fig.  1.  Hie  user  accesses  Tymshare,  Inc.  (with  suitable  pass- 
words, etc.)  and  asks  for  one  of  the  already  filed  sets  of  model  INPUT 
parameters  for  various  postures  and  situations  as  listed  in  the  "Catalog," 
given  later  herein,  or  creates  new  ones  from  scratch.  (Note,  the  "Catalog" 
files  are  "READ-ONLY"  so  they  cannot  be  inadvertantly  deleted.)  If  some 
of  the  parameters  need  to  be  changed,  the  user  employs  Tymshare' s EDITOR 
to  change  the  appropriate  numbers  on  given  lines  of  the  selected  "INPUT" 
file  (and  reflles  it  appropriately),  or  he  can  use  the  selected  set  without 
modification.  This  process  could  take  from  one  minute  to  one  hour  depend- 
ing on  the  number  of  parameters  to  be  changed  and  checked. 

The  user  then  identifies  the  desired  transfer  functions  via  a "CHOICES" 
file  wherein  each  desired  output/input  vibration  direction  is  listed.  This 
process  takes  from  1 to  10  minutes  depending  on  the  number  of  new  choices. 

The  user  then  activates  the  program  by  typing  "GO  BIODYN",  and  within 
seconds  the  process  described  above  is  executed  and  the  results  are  then 
typed  out  in  the  selected  format.  Selected  output  files,  so  identified 
in  the  "CHOICES"  file,  can  be  saved  for  later  plotting  or  recall.  From 
"GO"  to  printout  often  takes  less  time  than  it  did  to  read  this  paragraph. 
This  .computational  speed  makes  it  feasible  for  a user  to  sit  at  a terminal 
and  iteratively  alter  the  INPUT  and  CHOICES  files  to  iterate  on  the  charac- 
teristics of  a desired  biodynamic  transfer  function  or  mode* 

For  further  manipulation  of  the  resulting  transfer  functions,  such  as 
frequency  response  (Bode)  plots,  calculating  rms  output  in  response  to  an 
input  spectra,  etc.,  the  output  files  can  readily  he  accessed  by  the  user's 
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(or  other  Tymshare- resident)  analysis  programs.  A typical  set  is  the  STI 
Unified  Servo  Analysis  Method  (USAM)  programs  which  permit  a variety  of 
systems  computations  and  plots  for  mnltiloop  pilot-vehicle  systems  of 
this  type  (Ref.  3).  Although  such  extant  programs  are  not  an  intrinsic 
part  of  BIODYTT-78,  the  BIODYN  output  files  can  be  easily  transferred  to 
them  in  the  Tymshare  system. 

B.  MODEL  AID)  PARAMETERS 

1 . Model  Description 

Figure  2 (updated  from  Fig.  2 of  Ref.  l)  presents  the  biomechanical 
model  and  defines  many  of  the  necessary  parameters  that  describe  the  nominal 
(or  trim)  situation.  It  utilizes  a "homologous,"  or  life-like  representa- 
tion, of  the  major  body  segments  in  their  orientations;  simplified  to  a 
minimum  number  of  lumped  parameter  equivalents.  The  biomechanical  features 
include : 

• Semisupine  torso;  sliding  hip  plus  rocking  chest 
supported  on  a compliant  buttocks/seat. 

• Head  bobbing  on  an  articulated  neck  with  passive 
compliance. 

• Upper  arm  and  forearm  links  plus  grip- interface 
compliance,  driven  by  an  active  neuromuscular 
system. 

• Arm-rest  restraints  (optional). 

• Stick  "feel  system"  dynamics. 

This  simplified  model  was  derived  to  describe  the  motions  of  the  head 
and  arm  elements;  the  "pin- joint"  node  between  upper  and  lower  torso  segments 
Is  not  meant  to  represent  any  physical  feature. 

The  active  neuromuscular  system  shown  in  Fig.  2 is  a schematic  repre- 
sentation of  the  net  effect  of  the  complex  agonist/antagonist  muscle  pairs 
controlling  the  upper  arm.  A linearized  neuromuscular  system  representation 
is  shown  in  Fig.  3,  which  further  relates  this  to  the  effective  sensors  of 
muscle  length  and  force  as  well  as  propr.  oceptive  senses  from  the  stick 
grip  interface  which  completes  a loop  through  central  processing.  Unless 
the  neuromuscular  properties  are  being  investigated,  we  recommend  using 

k 


TR- 1037-2 


Figure  3. . Linearized  Neuromuscular  System 


the  given  values  which  are  representative  of  any  normal  person’s  arm-hand 
system,  and  generally  yield  reasonably  damped  neuromuscular  servo  proper- 
ties. Reference  k may  be  consulted  for  further  information  on  this  neuro- 
muscular model. 

2.  Parameter  Definitions 

Table  1 * presents  the  parameter  definitions  and  typical  values  used  in 
the  Ref.  1 validating  experiments.  This  is  an  updated  version  of  Table  A-i 
in  Ref.  1,  wherein  some  definitions  have  been  added  (e.g.,  stick  length, 
arm  rest  parameters)  and  a few  typographical  errors  have  been  corrected. 

The  equations-of-motion  for  BIODYN-78  have  been  written  in  a "second 
order  element"  matrix  form  using  generalized  (all  angle)  parameters.  When 
a particular  set  of  postural,  neuromuscular,  stick  and  display  dimensions 
are  input  as  a file,  then  the  program  internally  sets  up  the  general  equa- 
tions and  evaluates  the  linearized  coefficients  for  perturbations  about 
the  equilibrium  posture. 

C.  FILES 

1 . General 

Two  files  are  needed  to  execute  BIODYN.  The  INPUT  file  specifies  the 
80-odd  parameters  for  a given  model,  while  the  CHOICES  file  specifies  the 
transfer  functions  (output  variable s/forcing  function)  to-be  computed, 
printed,  arid  (optionally)  saved  bn  file. 

For  situations  involving  a seated  operator,  the  user  can  access  either 
of  the  two  model  cases  currently  resident  in  the  Tymshare  system  for  use 
in  BIODYN.  The  file  names  are  cataloged  below,  listings  later  herein. 

These  cases  were  used  to  validate  experimental  data  as  described  and  dis- 
cussed in  Ref.  1 (Appendix,  herein). 


*This  Table,  and  all  subsequent  Tables  are  located  at  the  end 
of  this  report. , 
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BIODYN-78  catalog  ' 


DESIGNATION  IN  RET.  1 (ALSO 
FILE  NAME  TABLE  l) 


COMMENTS 


B.19 


B.51 


"VIBRATION  EXPERIMENTS"  . Used  here  as  "test  case";  pilot 

in  standard  13  deg  seat,  hand 
on  center- stick 


"RIDE  QUALITY  EXPERIMENTS" 


Relaxed,  seated  operator,  fac- 
ing a display,  hands  in  lap 


These  files  are  used  in  various  ways,  depending  on  the  application. 

For  example; 

• Model  fitting  to  a given  set  of  data.  Select.,  a set 

of  transfer  functions  representing  the  measured  outputs/ 
inputs  (e.g.,  stick  or  head  motions  to  vertical  acceler- 
ation; . Create  a set  of  model  INPUT  files  bracketing 
the  likely  parameters,  at  .d/or  progressively  change  the 
INPUT  file  as  answers  are  inspected. 

• Choosing  measurement  points  for  a future  experiment. 

Create  one  model  INPUT  file  for  a representative  pos-. 
ture,  and  two  others  reflecting  the  likely  extremes., 

Create  a CHOICES  file  to  cover  each  likely  measurement 
point . 

Use  of  the  Tymshare ' s very  convenient  EDITOR  greatly  simplifies  the 
selection,  modification,  or  creation  of  the  model  parameter  files  or  output 
choices.  Sections  1-3  of  the  Tymshare  XEXEC  Manual  (Ref.  2)  should  be  read 
before  performing  these  functions. 

.2.  Model  INPUT  Files 

Table  2 documents  the  Catalog  "INPUT"  files  for  the  "Vibration  and  Ride 
Quality"  experiments  data  listed  in  Table  1.  Note  how  associated  groups 
of  parameters  in  the  vertical  arrangement  of  Table  • 1 are  changed  to  hori- 
zontal groups  of  parameters,  as  many  as  nine  numbers  in  a horizontal  line 
in  Table  2.  The  first  line  is  a "title"  line  (<  60  characters),  whereas 
the  contents  of  the  other  eleven  lines  are  indicated  by  the  separators  in 
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Table  1 . Finally,  in  Table  2 note  the  "free  form"  array  format,  -wherein  one 
or  more  blank  spaces  are  used  to  separate  the  numbers  in  a clear  manner. 

3.  Transfer  Function  CHOICES  File 


Table  3 summarizes  all  of  the  response  variables,  describing  the  various 
model  variables  of  Fig.  3,  as  well  as  element  perturbation  movements  arid  ' 
some  interface  forces  of  Fig.  2.  These  are  available  as  outputs  in  response 
to  the  inputs  (forcing  functions)  summarized  in  Table  k.  In  general,  the 
user  requests  transfer  functions  using  pairs  of  three-character  column 
codes  selected  from  Tables  3 and  k. 

The  transfer  functions  are  specified  in  a vertical  list,  putting  first, 
all  those  that  are  to  be  both  printed  out  and  saved  on  file  (user  supplies 
file  names  during  program  execution)  and  then,  below  a symbol  "XXX"  desig- 
nator, all  those  which  are  just  to  be  printed  out.  The  general  CHOICES 
file  format  is  thus: 


Select  from: 


Output/input: 


Table  Table 

3 

* * 

, AAA/BBB 

ccc/ddd 


To  print  and  save 
on  file 


xxx  Designator 

KRR/BBB  \ 

SSS/DDD  I 

. . \ To  print  only 


These  will  be  illustrated  in  the  Check  Case  example,  to  follow. 

The  response  variables  given  in  Table  3 for  shoulder  and  head  trans- 
latory  motion  are  defined  as  relative  to  the  platform.  The  inertial 
shoulder  and  head  responses  can  be  obtained  by  using  a set  of  "special" 
platform  deflection  forcing  functions  (denoted  by  BXP  or  BZP)  as  shown 
below:  • 
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CHOICES  FILE  ENTRY 


RESULT 


TRANSFER  FUNCTION 


inertial  shoulder  (vert . ) 
platform.  ~^ion  (vert.) 

inertial  shoulder  (horiz. ) 
platform  motion  (horiz.)  . 

inertial  head  (vert.) 
platform  .-notion  (vert.) 

IMPORTANT  NOTE:  ^ denptes  a typed  blank 

If.  Saved  Transfer  Fur ct ion  Files 

As  the  program  computes  each  numerator  which  is  to  be  saved  (those 
above  the  XXX  in  the  CHOICES  file),  it  stops  and  prompts  the  user  to  type 
in  a file  name  under  which  the  transfer  function  is  to  be  stored.  File 
naming  must  follow  Tymshare,  Inc.  conventions,  e.g.,  as  given  in  Section  3 
of  the  Tymshare  XEXEC  Manual  (Ref.  2),  In  general,  the  name  can  have  up 
to  six  characters,  followed  (optionally)  by  a period  and  up  to  three  more 
characters.  The  first  character  must  be  a letter;  the  others  can  be  any 
combination  of  numbers  and  letters;  A "temporary"  file  (automatically 
erased  at  the  end  of  a,  session)  is  designated  by  a .TMP  at  the  end  of  the 
file  name. 

For  those  transfer  functions  saved  on  file  the  following  format  is 
used;  FORTRAN  formats  are  shown  in  (parentheses),  for  those  who  may  wish 
to  read  the  files. 

NZ1,  NZ2,  NP1,  NP2,  K (5 X,  4(l3,lH,),  E12.6) 

1 Az-j 

V?z2 
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first  order  zeros  (one  per  line) 
(5X,  El 2. 6) 
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tfZl/BZP  fs 

ZP 

j&S/BXP  fs 

*P 

DZH/BZP  zh 


i 


1 


£z-|>  “Zj 

£z 2*  ^2 


second  order  zeros  (£2.,aj-,)  Per  line 

(aY . **13  £ r *7  fi?  £\ 


'/TP1 

l/^?2 


first  order  poles  (one  per  line) 
(5X,  £12.6) 


£pi>  “Pi 

C‘2'  ^2 


second  order  poles  ( Cpj^p)  Per  line 
(5X,  El 2.6,  1H,,  El 2 .6) 


•where : 

NZ1  = flush er  of  f:  'st-order  zeros 

NZ2  = Number  of  second-order  zero  pairs 

NP1  - Number  of  first-order  poles 

NP2  =»  Number  of  second-order  pole  pairs 

K = Gain  of  high-frequency  asymptote  (units  of  output/ 
' input  dimensions) 


An  example  is  shown  later. 

These  saved  files  may  be  used  later,  in  other  Tymshaxe-resident  programs, 
such  as  reformatting,  plotting,  incorporation  in  overall  system  models,  etc. 
A Bode  Plot  example  is  shown  later. 
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D.  COMPUTER  OPERATIONS  A!©  BIODYN  EXECUTION 
1 . Example  "Check  Case" 

Given  the  INPUT  and  CHOICES  files,  the  BIODYN  program  can  he  executed 
very  simply.  Figure  U is  an  example  of  actual  computer  operations  involved 
in,  first,  printing  out  files  to  be  used  (recommended,  practice  in  case  of 
errors),  and,  second,  actual  execution  of  BIGDYN-78.  This  check  case  uses 
file  B. 19  (the  Vibration  Feedthrough  experiments). 

The  underlined  items  in  Fig.  4 identify  those  typed  in  by  the  user 
(terminated  by  a "carriage  return")  — all  other  printing  was  done  by 
BIODYN  or  by  Tymshare's  operating  system  (called  XEXEC,  the  same  name  as 
Ref.  2) . XEXEC  prints  a prompt  (— ) at  the  beginning  of  a line  when  it  is 
ready  to  accept  a command  and  waits  there  for  a user  response.  The  circled 
numbers  on  the  left  in  Fig.  4 identify  items  discussed  below. 

0 User  issues  command  to  "TYPE"  out  "INPUT"  file  name  B.l9> 
followed  by  a carriage  return.  Then  XEXEC  lists  out  file 
contents. 

0 Same  for  "CHOICES"  file  (this  file  was  composed  earlier 
using  EDITOR,  Ref.  2. 

(3)  User  issues  "GO"  command  followed  by  file  name  "BIODYN" 
and  program  execution  begins. 

(X)  BIODYN  prompts  user  for  "INFJT"  file  name.  After  his 
response,  3.19>  program  sets  up  equations  of  motion. 

(5)  BIODYN  prompts  for  a "CHOICES"  file  name  and  U3er  types 
■ C.TMP.  Execution  continues  a3  program  prints  date  and 
tine  followed  by  the  case  title  (taken  from  the  first 
line  of  file  B.19  listed  above). 

0 The  denominator  i3  computed  and  printed.  The  first 
number  listed  is  the  gain  of  the  high  frequency  asynq)- 
tote. 

0 First  order  poles  are  listed,  up  to  fdur  per  line. 

0 Second  order  poles  are  listed  as:  (£,  u>,  £u>,  tu%/1  _ ^2) 

0 Gain  of  the  denominator  low-frequency  asymptote. 
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List  INPUT  file 


-TYPE  8, ,19. 
STIff  STICK 


1.  172 

. 29  . 

44 

.4124  44. 

1,417 

.-3*5 

. S 

.4152  145. 

.31 

2. 

1392*,  , 

.557345E-44  .. 

14. 

12 

92.  28 

29585. 27 

2.  1 

8*. 

4*. 

2.43123  144. 

5 . 

4939  4. 

*. 

.424437243 

.8  1 

6. 

i. 

.5 

,45555555  4. 

18. 

. 1 

5 .8 

1*. 

16.75  544. 

4.  34 

. . 

3 .439  5, 

. .44  1.  .4 

13. 

.3 

5. 

.75 

4.  . 

1 

4.  -24. 

-.25  54.  .45 

NOTE:  User  entries  are 
underlined 


-TYPE  C,T*9 
DTH/.DZ  P 

XXX 

C /AZP 

-GO  a:  coy*' 


ENTER  INPUT  PILE  NA*E:  B. 19  j 

ENTER  CHOICES  PICE  MANE:  C.TMP  ) 


ll-Jun-79  21:42 


CASE:  STirP  STICK 


DEMON  (MATOR  t 


List  CHOICES  file 


Activate  BIODYN 


Program  proapta 
for  file  names 


Case  Identifier 


.27577E-17 

16.593 

) 

( 58.788 

) l 97.343 

> > 

282.17 

> L 

357.45 
( .26496 

14.914 

8 

2.8917  , 

14.524 

)> 

r 

\ 

( .24252 

12.943 

2.6213 

12.675 

)) 

/ 

( .35427 

26.761 

9.4846  . 

25.425 

)) 

f 

( .15112 

29.611 

4.4747  , 

29.274 

( .85762 

32.953 

28.261 

16.947 

)> 

l 

< .19124 

8 

57.483 

t 

14. 993  . 

56.422 

)) 

/ 

TSASSFSR  FUNCTION  COMPONENTS: 

High  Frequency  gain 
{ Denominator ) 

first-order  poles 


< . 1187  3E»14> 


NUMERATOR:  DTN/DZP  riLE  MANE?  PTH.TMP 

NEW-  rue 


h 


12319E-18 

( .44444  ) 

( .44444 

) ( 16.594 

) ( 

18.21 

( 95.827  ) 

( 232.81 

) ( 498.67 

) 

((  .455*7  , 

14. 39* 

, 4.7392 

9.2528 

1) 

(I  .975541-41, 

16.243 

, 1.  5846  , 

16. 125 

)> 

((  .84471  , 

33.835 

, 28.445  . , 

18.321 

)) 

((  .23595  , 

51.815 

, 12.226 

54.  352 

>> 

((-.16361  , 

<-. 29394E*87> 

284.74 

,-75.464  . 

274.66 

)) 

I 

h 


Second-order  poles 

(?»  o>»  n>y,  _ ,i. 

Low  frequency  gain 
(Denominator) 


Program  prompt!  for 
file  name 

High  frequency  gain 
(Numerator) 

PI-  -order  zeros 


(Second-order  xeroe) 

Low  frequency  gain 
(Numerator) 


NUMERATOR  i C 
. 384274-18 


< 11.444  ) 

( 14. >37 

> 1 22.258 

) ( 

23.4 

(-38.117  ) 

< 114.44 

) 

(!  .19465  , 

12.  356 

9 

2.4(68  , 

12.147 

>> 

((  .65427  , 

17.818 

0 

15.222  , 

4.2626 

)) 

((  .14461  , 

28.944 

9 

4.2)81  , 

28.  591 

1 1 

((  : 974  41  , 

<-. 57)271. 14> 

147.45 

0 

143.34  . , 

26.428 

I) 

TRU  • 37.42 

ll-Jun-79 

21:43 

similar  to 
above 


Figure  Example  of  Computer  Operations  and 
BIODYN  Execution 
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(TO)  Biodyn  prompts  the  user  for  a file  name  to  store  the 
transfer  function.  User  responds  with  DTH.  TMP  and 
BIODYN  then  checks  the  user  file  directory  and  indi- 
cates that  this  is  either  an  OLD  FILE,  or  NEW  FILE  as 
here,  and  awaits  a user  response.  If  the  user  types 
a carriage  return  the  information  is  written  onto  the 
file,  destroying  the  previous  contents.  If  the  user 
enters  any  other  character  followed  by  a carriage  return, 
no  information  is  written  onto  the  given  file  name  — 
thus  sparing  the  user  from  destroying  possibly  valuable  ■ 
information.  Instead,  the  program,  prompts  RE-ENTER 
FILE  NAME.  After  user  responds,  the  program  indicates 
OID  FILE  or  NEW  FILE,  etc. 

m)  The  numerator  is  printed  out  following  the  same  format 
as  ©to®. 

(12)  The  next  numerator  in  the  CHOICES  file  is  printed  out. 

(13)  As  there  are  no  more  items  in  the  CHOICES  file,  the  pro- 
gram reports  the  number  of  TRUs  used.  ($.  13  each),  prints 
date  and  time  and  then  the  program  terminates  with  EXIT. 

The  dash  indicates  that  XEXEC  is  awaiting  a command.  At 
this  point  the  user  can  return  to  © or  (3)  , and  repeat 
this  entire  process  with  another  INPUT  and  CHOICES  file 
pair. 

Noce  that  each  transfer  function  is  to  be  interpreted  as 

<{n(s  + ^)}'{n(s2  + 2yoz*  + a|)| 

{<l(«  + ^)}*{ll(s2  + + 

where  II  denotes  a product  of  first-  or  second-order  root3,  and  < is  the 
so-called  "root  locus  gain"  (of  the  high  frequency  asymptote): 


High  Frequency  Gain  of  Numerator 
High  Frequency  Gain  of  Denominator 

and  the  various  first  and  second  order  poles  and  zero  are  indicated  in 
Fig.  4. 
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Transfer  Function  s - (s)  = 
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2.  Saved  File  Printout 


The  transfer  functions  that  are  saved  on  file  can  be  called  up  later 
by  other  programs  for  Bode  plotting  or  other  additional  calculations.  • 
Figure  5 is  an  example  printout  of  the  file  DTH.TMP  that  was  saved  on 
file  (Fig.  *0 . This  follows  the  format  given  earlier  in  the  subsection 
on  Saved  Transfer  Function  Files. 


— TYPE  DTH.TMP 

7,  5,  5,  6, 446729E-01 

.000000E+00 

.000000E +00 
* 165902E+02 
. 182  84  5E+02 
. 958270E-T02 
. 23281 5E-t-03 
. 498672Ef03 

. 455870E-r00,  . 103959E'*-02 
. 975536E-01 , .162027E+02 
. 840709E-r00 , .338348E+32 
. 235948E+00 , . 518149E-t-02 
-« 2636142+00 , .284736E+03 
. 165930E+02 
.587883E1-02 
. 973407E+02 
• 282173E+03 
. 357448E+03 

.264963E+00,  . 109137Ef02 
. 202525E+00,  . 1294  29E+02 
.354269E*00,  . 267610Et02 
. 151120E+00,  . 296106E+02 

• 857619E+00 , * 329528E-I-02 

• 191240E-*-00 , . 574830E-f02 


Figure  5.  Example  Printout  of  a Saved  File 
(User  entries  are  underlined) 


3.  Error  Massages 


There  are  two  types  of  user  errors  that  can  trigger  various  messages  — 
spelling  errors  in  the  CHOICES  file  and  numerical  errors  in  the  INPUT  file. 

a.  Spelling  errors.  In  the  CHOICES  file  a misspelled  response  variable 
or  forcing  function  cannot  be  detected  by  the  program  unless  it  results 
in  a three-character  code  that  is  not  in  Table  3 or 

Thus,  for  example,  a transfer  function  request  of  DHZ/BZP  cannot  pro- 
duce a numerator  since  DHZ  is  not  a response  variable  (Table  3).  In  this 
case,  the  following  message  will  be  printed  where  the  numerator  would 
have  been  printed 


DHZ  NOT  FOUND'  DHZ  BZP  SWAP 


CHOICES  file  requests 

The  same  type  of  message  results  for  a not-found  forcing  function,  i.e., 
the  "not  found"  characters  are  always  printed  first  followed  by  both  por- 
tions of  the  original  request. 

b.  Numerical.  A number  of  error  messages  can  result  from  grossly  inappro- 
priate numerical  selections.  The  equations  of  motion  have  been  written  in 
a manner  such  that  most  masses  and  many  other  elements  can  be  set  to  zero, 
e.g.,  the  many  zeroed  elements  for  the  Ride  Qpality  Experiments  (Table  i). 
However,  if  the  user  also  zeros  up  Kin  or  or  Kp  then  the  error  message 

NULL  MATRIX 

will  appear  in  the  output  in  place  of  the  DENOMINATOR.  For  other  numerical 
selections  this  message  might  appear  among  some  of  the  numerators.  In 
general,  the  user  should  examine  the  INPUT  data  to  check  for  typing  errors, 
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e.g.,  a massless  unrestrained  head  such  that  Aejj  is  indeterminate  (NULL 
MATRIX  will  result). 

A number  of  other  error  messages  of  varying  degrees’  of  severity  and 
user  response  can  occur.  OVERFLOWS  and  UNDERFLOWS  are  printed  by  the 
operating  system.  If  only  one  of  each  of  these  occurs,  they  probably 
occur  in  intermediate  calculations  and  can  often  be  ignored  if  an  other- 
wise normal  printout  (e.g.,  Fig.  1)  follows,  since  roots  are  not  accepted 
until  they  satisfy  a convergence  criteria. 

Some  of  these  problems  can  be  solved  by  judicious  parameter  selection, 
viz.,  in  the  Vibration  Experiments  (Table  l)  the  neck  mass  and  inertia  were 
set  to  zero.  The  latter  caused  no  problems  in  the  equations  of  motion, 
whereas  very  small  values  would  produce  an  additional  quadratic  mode  at 
very  high  frequencies  that  would  appear  both  in  the  denominator  and  all 
numerators.  These  high  frequency  "nuisance  modes"  are  not  important  to 
any  system  responses,  yet  keeping  them  in  forces  the  root-finding  routine 
to  hunt  them  down  with  attendant  OVERFLOW  messages  if  they  are  greater  than 
about  1000  rad/sec.  Similar  problems  can  result  for  very  small  values  of 
some  dampers,  again  pushing  roots  to  very  high  frequencies  such  that  the 
following  message  might  appear: 

EXCESS  OR  INCORRECT  ROOT  POUND 

Thus  it  is  better  to  avoid  using  very  small  values  of  masses  and  dampers 
for  those  elements  not  important  in  the  user's  application.  The  authors 
urge  the  potential  user  to  study  the  differences  between  the  parameter 
sets  (Table  l)  used  in  the  Vibration  Experiments  (17th  order  denominator) 
and  the  Ride  Quality  Experiments  (9th  order  denominator).  The  variety  of 
elements  "zeroed  up"  in  each  study  helped  remove  "nuisance  modes"  such  that 
there  were  no  error  messages. 

Finally,  to  prevent  "runaways"  (rare  but  potentially  costly),  the  follow- 
ing message  may  occur: 

. 50  ITERATIONS  WITHOUT  CONVERGENCE,  CONTINUE?  (Y  OR  N): 
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The  user  responds  with  either  Y or  N,  followed  by  a carriage  return. 

The  first  response  continues  the  root  search  for  up  to  50  more  iterations 
while  the  latter  accepts  the  roots  found  so  far.  If  no  OVERFLOWS  have 
occurred  the  authors  recommend  responding  with  Y since  the  user's  para- 
meter set  may  have  an  unusual  Root  configuration  and  allowing  the  program 
to  continue  may  produce  convergence  for  all  roots  before  this  message 
occurs  again.  Responding  N followed  by  carriage  return  will  stop  the 
search  and  report  out  all  roots  found  so  far,  but  not  all  that  may  be  pre- 
sent. BIODYN  will  continue  normally  with  the  next  item  in  the  CHOICES  file. 
A special  investigation  found  that  each  root  found  ■por  the  Vibration  Experi- 
ments or  the  Ride  Quality  Experiments  required  much  less  than  50  iterations. 

4.  Example  Bode  Plot  of  a Saved  File 

Figure  6 illustrates  the  use  of  STI’s  TYMSHARE  resident  Unified  Servo 
Analysis  Method  (USAM)  program  (Ref.  5)  to  access  the  transfer  function 
stored  on  file  name  DTH.TMP  (Fig.  4)  and  subsequent  terminal  Bode  plot. 

While  this  USAM  program  is  not  a part  of  BIODYN,  it  is  readily  available 
for  this  coarse  printer  plot,  as  well  as  better  resolution  ZETA  graphics 
Ref.  3 documents. 

J.  U3e  in  PIVTB  Program 

The  performance  computation  program  "PIVIB"  described  in  Ref.  5 requires 
various  biodynamic  feedthrough  transfer  functions  to  be  supplied,  and 
BIODYN-78  is  appropriate  for  this  purpose.  BIODYN  computes  the  X,  Z,  and 
9 outputs  required  in  pole,  zero,  and  gain  form  (designated  POLE0  in  PIVIB). 

For  example,  for  the  feedthrough  to  the  control  stick  from  vertical  or 
fore-aft  vibration  the  CHOICES  file  entries  would  be  as  follows : 
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♦Reference  5,  Paragraph  7*2.3 

Similar  procedures  are  used  for  any  other  biodynamic  transfer  function 
required  in  the  BIOTR  Submodule  in  PIVTB.  While  BlODYN  uses  the  international 
standard  units  (SI  system)  throughout,  PIVTB  uses  a variety  of  english  units, 
and  so  gains  and  dimensions  must  be  scaled  accordingly. 
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User  entries  are  underlined: 

Selects  file-name  • 
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Figure  6.  Typical  USAM2  Terminal  Bode  Plot  for 

Check  Case 
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TABLE  1 . PARAh&JCER  DEFINITIONS  AND  TYPICAL  VALUES 


PARA-  FORTRAN 
PIETER  NAME 


deg  8, 


kg  H?  XM2 

a Lg  B2 

a/a  Uj  Lg  R2 

kg*ffl?  Ig  Z2 

deg  8e  TE 


0 kg  Mg  XMS 

6 N/m/s  B,  BS 

0.  N/a  Kj  XKS 

'558  m/N  Cx  Cl 


0.0 

0.’6i 
13900.  N/a 


Bx/Kj  TIF 
. XLC 


14.  kg  Mg  XMB 

315-  n/m/s  Bj  , BB 

1*9360.  N/a  Kg  XKB 


•Source  of  data, In  Appendix  A 
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File  Name 


Caae  Identifier 


' Upper  arm  caaa 

Upper  am  length 
(*  * Lb> 

Upper  am  c.g. 
(fractional  distance) 

Upper  am  inertia 

Upper  am  angle 


lower  am  mass 

lower  am  length  (=  Le  + Ig) 

lower  am  c.g. 

(fractional  distance) 

lower  am  inertia 

Elbow  angle 


Stick  mass 
(referred  to  grip) 

Stick  damper 

Stick  gradient 

Grip  interface 
compliance,  K^* 

Grip  Interface  tine 
constant 

Stick  length  (piTOt  to 
center  of  grip) 

Scale  factor  (multiplies 
stick  displacement) 

Stick  compliance  parameter 
(CKS  » 0 means  that  Kjg  » w 


Lower  body  nass  (Hips) 

lower  body/seat  cushion 
daoper 

Lower  body/seat  cushion 
stiffness 


Control 
Stick  and 
Grip  Inter- 
face 


Lower  Bodv 


TYPICAL  VALUES 


vi3  ration 
EXPERIMKTS* 

rite  quality 

EXPERIMENTS* 

UN  TPS 

18. 

18. 

kg 

0.15 

0.3 

m 

0.8 

0.8 

kg-m2 

10. 

-13- 

deg 

16.75 

Uo. 

N-m 

rad/s 

500. 

500. 

N-m 

rad 

U.34 

3-1 

kg 

0.3 

0.15 

a 

0.039 

0.0303 

kg-mS 

5- 

-5- 

deg 

0.0 

0.126 

N*m 

rad/s 

1. 

15-  . 

0.0 

0.025 

PARA- 

METER 

FORTRAN 

NAME 

DEFINITION 

XMT 

Torso  mass 

Lr 

XLT 

Tbiso  c.g.  to  hips 
pivot  length 

zrr 

Torso  inertia 

B 

TET 

Angle  of  XLT  from 
vertical 

Torso 

BTB 

Torso  damping 

, 

KtB 

XKTB 

Torso  stiffness 

XMi 

Head  mass 

ITN 

XLTN 

Torso  c.g.  to  neck/ 
torso  pivot 

lH 

ZIH 

Head  inertia 

eTN 

raw 

Angle  of  XLTN 

®HH 

BHN 

i 

Kead/neck  damper 

% 

| 

XXHN 

Head/neck  stiffness 

Head/Neck 

% 

XLH 

1 

Head  c.g.  to 'head/ 
neck  pivot  length 

®H 

TH 

Angle  of  XIH 

cH 

CH 

• 

Head/heck  compliance 
parameter  (CH  =*  0.  locks 
had  on  neck) 

9L 

1 

THL 

Seat  (hips)  tilt  angle 

■ 

L. 

i 

XLS 

Tbrao  c.g.  to  shoulder 
length 

Torso 

Js 

THS 

Angle  of  XLS 

Jb 

VD 

Viewing , distance 

Display 

bar 

BAR'  ■ 

Arm  rest  damper  (normal)  . 

^AR 

XKAR 

Arm  rest  stiffness 
(normal) 

Sp 

BAT 

Arm  rest  damper 
(tangential) 

XKAT 

Arm  rest  stiffness 
(tangential) 

Arm 

Rest' 

XLER 

Elbow  to  arm  rest 
distance 

ARMS 

■ 

Fraction  of  arm  weight 
on  arm  rest  j 

I 


TABLE  1.  (Concluded) 


* FILE 

— 

typical  va lues 

PAHA- 

METER 

FORTRAN 

NAME 

DEFINITION 

. lh r« 

NO. 

VI3HAT  ion 
.EXPERIMENTS* 

RIDS  QUALITY 
EX?ERIME!CS* 

UNITS 

12 

0.0 

1 .24 

kg 

% 

XMK 

Neck  mass 

0.1 

.12444 

tt 

% 

na 

Neck  length 

0.0 

0.0024 

kg-n2 

% 

zvn 

Neck  inertia 

. 

-20. 

-50. 

deg 

* 

THN 

Angle  of  XLN 

Neck/Torso 

0.25 

0.0 

N-m 

BPT 

Neck/torso  damper 

rad/? 

50. 

50. 

KffT 

XKWT 

Neck/ torso  stiffness 

0.05 

.06222 

XLNl 

Neck  c.g.  to  neck/ tor so 

pivot 

-50. 

-50. 

■ 

deg 

9v 

THV 

Angle  of  line  of  sight 

Display 

TABLE  2.  LISTING  OF  BIQDYN-78  CATALOG  FILES 


VIBRATION  EXPERIMENTS 
(file  name  B.19) 


STIFF  STICK 

1.372  .29  .44  .0120  40. 

1.017  .305  .5  .0152  145. 

.31  2.  13900.  . 557305E-04  .01  .61  13900.  1. 

14.  1290.08  29585.27 

2.  1.  80.  40.  2.43123  104.39969  1.  71519.  0. 

5.  .0909  0.  0.  .020437203  .09090909  .08988764 

.8  16.  0.  .5  .05555555  0.  0.  1. 

18.  .15  .8  10.  16.75  500. 

4.34  .3  .039  5.  .00  1.  .0  70.  0. 

13.  .3  5.  .75 

0.  .1  0.  -20.  .25  50.  .05  -30. 


RIDE  QUALITY  EXPERIMENTS 
(file  name  B.51 ) 

HFGn.NH'CK  ••  rOf-'SO  ONLY 

>.ooo  .44  .0000  40. 

o.  . u .oooo  145. 

.00  0.  13900.  .557805E-04  .00  .61  13900.  1. 
14.  31V, . 49360. 

’•  0.  30.  40.  0.  0.1.  119322.  0. 

0. 

0.  ‘ , . 

i a.  .30  .8  -in.  40.  500. 

3.1  . r,  .0303  -5.  .126  15.0  .025  60.  1. 

13.  .3  -5.  .6853 

l .24  • 12444.  .0024  -30,  .00  50.  .06222  -30, 
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TABLE  5.  RESPONSE  VARIABLES 


TABLE  3.  (CONCLUDED) 
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•Note:  Snail  "(<"  represents  a tvped-blank  space 


t These  are  used  to  get  transfer  functions  of  Platform  notion 


